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ABSTRACT 

The  U-  S.  Naval  Civil  Engineering  Laboratory  is  studying  the  design  of 
various  entranceways  for  protective  structures.  The  particular  portion  of 
the  work  reported  herein  is  con  erned  with  air  ducts  in  concrete.  One  of  tl 
essential  parameters  for  the  design  optimization  of  such  ducti  is  tiie  neutri 
energy  spectra  as  a  function  of  position  in  the  duct.  From  this  parameter  r 
be  determined  the  type  and  quantity  of  shielding  materials  necessary  to  ach' 
a  specified  protection  factor- 

•  • 

In  order  to  achieve  a  uniform  system  of  neutron  dosimetry,  it  was  decir 
to  use  the  activation  of  bare  and  cadmium  covered  foils,  foil  sandwiches,  ai 
threshold  foil  detectors  to  measure  flux  spectra.  This  application  represei 
one  of  the  first  times  that  foils  were  used  to  experimentally  determine  dif¬ 
ferent  I a  1  flux  over  the  energy  range  from  0.001  ev  to  15  Mev.  It  is  believ* 
to  be  the  first  application  to  protective  structure  shielding  experiments. 

Preliminary  to  the  experimental  determination  of  neutron  spectra  by  fo 
activation,  it  is  essential  that  the  responses  of  the  selected  foils  be  mea: 
in  a  medium  of  knoym  moderating  and  scattering  properties.  A  hydrogenous  n* 
such  as  paraffin,  meets  this  requirement  and  is  experimentally  simple.  A  s< 
of  nominally  14-Mev  neutrons  was  directed  at  successively  increasing  thickn< 
of  paraffin  made  from  1-inch  slabs  12x24  inches.  Selected  foils  were  placet 
the  opposite  side  of  the  paraffin  and  their  response  to  the  neutron  irradia' 
determined.  From  this  study  technology  is  developed  by  which  the  observed  i 
ponse  of  foils  Irradiated  in  complex  shields  may  be  related  to  neutron  ener 
flux.  The  report  of  this  technology  will  be  made  at  the  time  neutron  erarg 
spectra  for  air  ducts  in  concrete  is  reported.  Herein  is  the  report  of  the 
preliminary  experiment  where  a  study  was  made  of  neutrons  streaming  through 
hydrogenous  media.  The  normalized  activation  responses  of  selected  foils  a 
given  as  a  function  of  paraffin  thickness. 


Also  described  fs  the  sequence  of  computer  programs  which  process  and 
analyze  the  radiometric  data  of  such  foils.  These  programs  were  specifically 
designed  to  handle  the  foil  data  associated  with  the  experimental  determination 
of  neutron  streaming  through  ducts;  however,  they  are  suitable  for  handling 
other  radiometric  analysis  problems.  By  their  use,  radiometric  analysis  data 
may  be  rapidly  processed  with  reasonable  assurance  of  freedom  from  numerical 
error.  Maximum  utility  of  these  programs  is  obtained  when  using  automatic 
counting  that  punches  data  directly  on  cards  usable  as  the  computer  input. 

This  work  was  sponsored  by  the  Defense  Atomic  Support  Agency. 


Qualified  requesters  may  obtain  copies  of  this  report  from  DDC. 
Release  to  the  Clearinghouse  is  Authorized. 


INTRODUCTION 


The  application  of  foil  activation  for  neutron  flux  measurement  in  j 
the  energy  range  above  1  Mev  and  below  0.5  ev  has  been  described  by  Hughes. 
Above  1  Mev  threshold  reactions  are  used  and  below  0.5  ev  bare  and  cadmium 
responses  are  used.  From  the  many  foil  techniques  available  for  the  measure¬ 
ment  of  neutrons  with  energy  between  0.5  ev  and  1  Mev,  the  application  of  foi 
sandwiches  was  selected.  This  application  has  been  described  by  the  a 
after  a  method  of  Stehn^  and  the  Staff  of  Argonne  National  Laboratory. 


The  techniques  for  calculating  neutron  energy  spectra  from  the  activation 
response  of  several  foils  has  been  previously  discussed.  Essential  to  such 
calculations  is  a  knowledge  of  the  various  responses  in  a  medium  of  known 
moderating  and  scattering  properties.  For  this  purpose  a  neutron  generator 
employing  the  d,t  reaction  with  an  acceleration  of  150  kev  was  used  as  a  sourc 
of  nominally  14-Mev  neutrons.  This  source  was  directed  at  a  praffin  moderator 
so  that  the  beam  of  deuterons  was  perpendicular  to  the  exterior  surface  of  the 
moderator.  The  response  of  selected  foils  was  observed  as  a  function  of 
paraffin  thickness  in  such  a  manner  that  no  reflection  was  present.  Each  ir¬ 
radiation  was  normalized,  maximizing  the  deuteron  beam  current  by  adjusting 
the  focus,  extraction,  and  accelerating  voltage,  and  observing  the  counting 
rate  ofaslightly  moderated  BF,  detector  used  as  a  flux  monitor.  The  high 
voltage  was  adjusted  last  and^assured  that  the  incident  neutron  energy  re¬ 
mained  constant  throughout  the  experiment.  This  is  possible  since  the  d,t 


reaction  peaks  at  approximately  150  kev  acceleration  at  which  value  the 
ejected  neutrons  have  a  most  probable  energy  of  14  Mev.' 


The  activity  of  each  foil  was  observed  with  an  internal  gas  flow  pro¬ 
portional  detector  with  a  "micromil"  window  and  using  "P"  gas  (90  percent 
argon,  10  percent  methane}.  The  detector  was  incorporated  into  an  auto¬ 
matic  counting  system  that  directly  punched  cards.  Data  was  identified  by 
separate  problem,  sample,  counting  system,  counting  geometry,  and  observa¬ 
tion  numbers.  All  foils  of  the  same  type  were  of  the  same  size  and  approx imat 
mass  and  were  irradiated  for  the  same  time.  Data  was  corrected  for  detector 
resolution  loss  and  decay  after  irradiation.  Data  was  normalized  in  such  a 
manner  that  knowledge  of  the  distance  from  foil  to  the  tritium  target, 
counting  efficiency,  foil  cross  section,  and  irradiation  time  was  not  part 
of  the  analysis. 


EXPERIMENT 


The  experiment  consisted  of  irradiating  certain  selected  foils,  listed 
in  Table  I,  with  neutrons.  Various  thicknesses  of  paraffin  were  placed  be¬ 
tween  the  neutron  source  and  the  foils  so  as  to  moderate  their  initial  nominal 
14-Mev  energy  and  to  cause  the  neutrons  to  take  on  differing  energy  distributi 
No  paraffin  was  placed  behind  the  foils  so  as  to  minimize  response  to  neutron 
backscatter.  For  each  specific  foil  material  the  same  foil  dimensions  were 
used  throughout  the  series  of  experiments;  however,  variations  occurred  betwee 


different  foil  materials.  The  neutrons  were  produced  by  accelerating  deuterium 
to  approximately  150  kev  before  impinging  on  a  tritium  target.  The  accelerator 
operating  conditions  were  optimized  throughout  the  experiment  by  adjusting  beam 
focus,  extraction  and  acceleration  voltage  so  as  to  maximize  the  beam  current. 

It  should  be  noted  that  the  d(t,n)  cross  section  is  near  its  maximum  value  at 
150  kev  energy.  At  this  energy  of  bombardment,  the  neutrons  have  a  most  probable 
energy  of  about  14  Mev.  All  foils  for  threshold  detectors  were  irradiated  30 
minutes.  The  aluminum-manganese  foil  sandwiches  were  irradiated  for  10  minutes 
followed  by  a  10-minute  period  for  manipulation  and  counting.  This  in  turn 
was  followed  by  a  15-minute  irradiation  for  all  cooper  foil  sandwiches.  The 
bare  and  cadmium  covered  indium  foil  combinations  were  irradiated  at  the  same 
time;  they  were  placed  for  irradiation  to  begin  with  the  10-minute  irradiation 
of  the  aluminum-manganese  sandwich  and  were  removed  after  the  15-minute  copper 
sandwich  irradiation.  Since  the  irradiation  time  of  each  specific  foil  was 
the  same,  no  corrections  were  made  for  activity  buildup  during  irradiation. 

A  slightly  moderated  BF,  detector  was  placed  in  a  fixed  position  near  the 
generator.  The  count  or  this  detector  integrated  over  the  irradiation  period 
was  used  as  a  relative  measure  of  neutron  flux  which  In  turn  was  used  to 
normalize  all  observed  foil  activities.  No  corrections  were  made  for  very 
slight  differences  in  mass  between  different  foils  of  the  same  material  or 
for  counter  efficiency.  However,  inverse  square  corrections  were  made  for 
the  varying  distance  between  the  foils  and  the  tritium  target.  All  observed 
foil  activities  were  corrected  for  decay  between  the  end  of  irradiation 
and  the  time  of  observation.  Particular  attention  was  directed  toward  preci¬ 
sion.  Care  was  taken  to  insure  the  purity  and  cleanliness  of  foils. 

Foil  activities  were  observed  with  an  internal  gas  flow  proportional 
detector  having  a  micromil  window.  This  detector  was  placed  in  an  automatic 
sample  changer.  The  counting  gas  was  90  percent  argon  and  10  percent  methane, 
and  counts  were  observed  at  the  center  of  the  beta  plateau.  A  RaDEF  source 
electro  deposited  on  a  platinum  disk  supplied  by  the  National  Bureau  of 
Standards  was  used  as  a  long  lived  standard  to  assist  in  obtaining  uniform 
counter  response  during  the  period  of  the  experiments.  The  output  of  the 
counter  was  the  input  to  a  data  converter  which  punched  on  IBM  cards  the 
counting  data.  Other  information  necessary  to  the  experiments  and  analysis 
of  data  was  manually  punched  on  these  and  other  cards.  A  complete  description 
of  the  data  format  on  these  cards  and  the  information  that  they  require  will 
be  found  in  che  following  section:  Data  Reduction.  All  counting  data  was 
analyzed  by  the  methods  described  therein,  which  included  statistical  analysis 
of  variance. 

DATA  REDUCTION 

Data  experimentally  obtained  from  automatic  counting  systems  during 
experiments  that  measure  neutron  streaming  through  ducts  may  best  be  reduced 
by  computer  techniques.  Computer  programs  are  also  suitable  for  handling 
other  types  of  radiometric  analysis.  The  processing  of  data  by  a  digital 
computer  insures  accurate  reduction  and  rapid  handling,  thus  affording  im¬ 
mediate  review  between  experiments.  The  programs  about  to  be  presented  were 


used  to  reduce  the  data  from  the  experiment  just  discussed.  However,  to 
Illustrate  the  broader  applications  of  these  programs,  numerical  examples 
are  included  from  reduction  of  other  counting  data. 

All  output  data  from  these  programs  may  be  used  as  data  with  which  to 
calculate  neutron  flux.  The  calculat IonsRof  this  program  are  performed 
basically  as  has  already  been  delineated.  Only  departures  from  previous 
work  and  new  Information  is  presented  herein.  The  programs  described  are: 
pre-  and  post-data  compiler,  PDC  1,  II;  three  data  organizations,  DATORG  I, 

II,  III;  resolution,  RESOL;  two  decisions,  OCISON  I,  II;  Chi  Square,  CHISQ; 
three  analyses  of  variance,  VAR  I,  Prep  VAR,  VAR  II;  three  decay  analyses, 
DECAY  I,  II,  III;  dead-time  determination,  DEAD;  and  backing  thickness 
determination,  BACKTH.  The  latter  two  programs  are  not  part  of  the  routine 
analysis,  but  are  included  for  completeness. 

The  program  sequence  Is  shown  by  the  block  diagram  of  Figure  1.  All 
programs  were  written  in  Fortran  I  and  were  designed  to  be  compiled  with 
the  PDQ  processor,  using  the  IBM  1620  computer.  This  computer  was  selected 
since  it  is  available  near  the  experimental  facility.  It  is,  however,  possible 
to  modify  the  program  for  use  on  larger  computers  (e.g.,  7090  or  709^  as  a 
chain  job  with  a  control  program).  This  would  be  accomplished  by  changing 
a  few  input  and  output  statements. 

PDC  takes  the  counting  data  as  punched  by  automatic  counting  equipment, 
whose  format  is  shown  in  Figure  2,  preceded  by  a  deck  of  program  definition 
cards.  The  format  of  the  problem  definition  cards  is  shown  in  Figure  3-  The 
following  programs  are  designed  to  be  run  in  sequence  according  to  on-line 
comments.  For  specialized  cases  of  radiometric  data  certain  of  the  programs 
could  be  omitted  from  the  running  sequence;  however,  the  experimentor  is 
cautioned  that  such  omission  could  result  in  a  data  format  that  the  computer 
does  not  properly  recognize. 

PDC  organizes  all  data  by  problem  number  and  sample  number  and  into  a 
format  most  suitable  for  further  reduction.  Its  output  is  the  input  to  DATORG 
which  sorts  the  data  according  to  counting  systems,  and  geometry.  These 
various  sortings  allow  the  intermixing  of  data,  thus  affording  flexibility  of 
counting  operations.  Its  output  is  the  Input  to  RESOLN  which  corrects  observed 
counts  for  detector  dead  time.  Data  is  then  fed  to  DC  I  SON  which  corrects  obser 
counts  for  decay  during  the  counting  period.  DCISON  also  tests  to  determine 
if  decay  has  occurred  between  two  successive  observations  on  the  same  sample. 

If  so,  each  observation  is  called  "number  one",  but  if  not,  the  observation 
number  is  advanced  sequentially,  regardless  of  the  observation  number  punched 
on  the  card.  The  output  from  this  program  is  the  input  either  to  DECAY  or  to 
CHISQ,  according  to  a  predetermined  code  number. 
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DECAY  corrects  for  decay  between  count  observations  and  adjusts  these 
observations  to  that  estimated  for  a  previous  time.  CHISQ  is  customarily 
used  only  when  a  new  system  is  put  into  operation  or  its  performance  is 
questioned.  It  determines  if  a  series  of  observations  is  from  the  Poisson 
distribution  associated  with  the  characteristic  randomness  of  radioactive 
disintegration.  Correct  performance  is  exhibited  only  if  such  a  distribution 
exists.  The  output  of  DECAY  is  the  input  to  VAR  which  averages  the  counts, 
corrected  to  a  particular  time,  that  are  associated  with  like  samples.  The 
error  associated  with  this  average  is  also  estimated  along  with  the  number 
of  additional  samples  necessary  for  a  smaller  error.  Not  directly  associated 
with  the  data  processing,  but  included  for  completeness,  are  two  programs, 

DEAD  and  BACKTH.  The  former  calculates  the  counter  dead  time  used  in  RESOL, 
and  the  latter  calculates  the  thickness  of  sample  backing  material  by  three 
different  method's:  first,  from  data  given  in  milligrams  per  square  centimeter; 
second,  from  data  given  in  mils  along  with  the  material  density;  third,  from 
the  ratio  of  the  observed  counts  of  another  source  with  and  without  the  backing 
material . 

The  input  and  output  formats,  flow  charts,  and  program  listings  will  be 
presented  later  herein  when  each  of  the  several  programs  is  considered  in 
more  detai 1 . 

Program  Description  -  Pre-Data  Compiler 

PDC  sorts  the  counting  data  problem  number,  the  major  data  breakdown,  and 
transforms  the  raw  data,  punched  by  the  automatic  counting  equipment,  into  a 
format  more  suitable  for  further  reduction.  The  first  cards  read  are  all  of 
the  constants  and  problem  identifying  numbers  required  for  further  data  re¬ 
duction,  such  as  sample  numbers  that  should  be  compared,  background  data,  decay 
constants,  number  of  decay  components,  and  special  code  numbers.  Next,  the 
data  furnished  by  the  automatic  counting  equipment  is  read.  It  should  be  noted 
that  the  sample  counting  time  does  not  have  the  decimal  point  included  on  the 
card  (a  requirement  of  the  counting  equipment).  Thus,  this  4ntry  is  immediately 
transferred  from  a  fixed  point  number  (no  decimal  point)  to  a  floating  point 
number  (with  decimal  point).  With  the  particular  counting  equipment  utilized, 
the  omitted  decimal  point  on  the  input  data  card  necessitates  dividing  the 
punched  number  by  100  in  order  to  obtain  the  counting  time  in  minutes.  After 
preliminary  sorting  of  data  into  separate  problems,  the  input  constants  on 
cards  are  punched,  followed  by  the  data  cards  as  they  are  processed.  These 
comprise  the  input  to  post-data  compiler,  which  organizes  the  problem  by 
sample  number  and  similarly  punches  the  data  in  the  same  format.  Table  II 
lists  the  source  program;  the  input  and  output  data  for  a  sample  problem 
are  given  in  Tables  1 L I a  and  I  I  lb,  respectively. 
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Data  Organization 

DATORG  consists  of  a  sequence  of  programs  beginning  with  the  output  of 
PDC.  It  arranges  the  data  In  ascending  order  with  respect  to  system  number, 
and  geometry  number,  associating  with  each  group  the  appropriate  Input  constan' 
The  Input  data  Is  punched  In  time  sequence  so  that  each  field  will  be 
organized  In  time.  If  the  data  Is  already  organized,  this  sequence  of 
programs  may  be  bypassed  and  the  output  of  PDC  used  directly  as  Input  to 
RESOL.  Once  the  data  has  been  stored,  a  loop  In  the  program  performs  the 
task  of  making  the  necessary  changes  to  the  sequence  of  data  cards  for  their 
organization.  After  the  data  Is  sequenced  appropriately  In  the  storage  array, 
duplicate  cards  are  punched  in  an  Identical  sequence.  At  this  point,  the 
next  set  of  raw  data  is  loaded,  to  be  arranged  In  similar  manner.  Tables  VI, 
VII and  X  list  the  source  program,  and  the  output  listing  of  a  sample 
problem  Is  shown  in  Tables  VII,  IX,  XI,  for  DATORG  I,  II,  arid  III. 

Resolution 

RESOLN  corrects  the  observed  counts  for  detector  dead  time  loss.  It 
should  be  noted  that  if  the  counting  rate  of  the  data  is  sufficiently  low 
so  as  not  to  exhibit  the  effects  of  resolution  loss,  this  program  may  be 
bypassed,  and  the  data  used  as  Input  directly  to  DCISON.  RESOL  obtains  the 
dead  time  from  the  data  cards  associated  with  the  corresponding  system.  After 
processing  data  for  a  particular  system,  a  new  system  number  and  value  of 
resolution  time  is  associated  from  the  input  with  the  counting  data.  The 
analysis  is  then  repeated.  Table  XII  is  the  source  program  listing  and 
Table  XIII  lists  the  output  for  a  sample  problem. 

Determining  Dead  Time 

Dead  time  is  calculated  from  a  series  of  count  rate  observations  made 
on  a  highly  active,  rapidly  decaying  radioactive  sample.  The  analysis 
utilizes  two  data  sets,  N,  the  true  count  rate  obtained  by  extrapolating 
low  count  rate  observations  with  a  known  half-life,  and  M‘,  the  observed 
high  count  rate  observations.  For  every  point  M‘  there  must  be  an  associated 
point  N.  A  curve  of  the  form  for  dead  time  corrections 

N  l-M'T 

where  T  is  dead  time,  is  then  fitted  to  the  data  by  the  least-squares  method. 
The  program  is  limited  to  100  data  points  for  each  of  the  M'  and  N  sets. 

Input  data  is  obtained  from  DCISON  I,  bypassing  RESOLN.  Table  XIV  is  a 
listing  of  the  source  program,  and  the  output  data  is  listed  in  Table  XV. 
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Dec  is  ion 


DCISON  determines  whether  or  not  radioactive  decay  is  significant.  It 
also  supplies  background  count  data  and  the  decay  constant  of  principal  intersect 
There  are  two  parts  of  this  program:  One,  DCISON  I,  is  used  when  the  sample 
observed  contains  only  one  radioisotope  (i.e.  decays  with  a  constant  half-life) 
and  the  other,  DCISON  II,  is  used  with  a  mixture  of  radioisotopes.  Input  data  is 
obtained  from  RESOLN  which  senses  the  correct  program  and  prints  on-line 
instructions  to  the  computer  operator.  If  no  background  count  information  is 
given  along  with  the  counting  data,  the  program  searches  for  this  information 
in  the  experimental  constants  section  of  the  input  data.  It  then  determines  if 
a  decay  constant  is  given  by  the  counting  data  cards.  If  not,  an  effective 
decay  constant  is  estimated  by  the  relation 


where  A.  and  A.  +  j  are  two  observations  of  count  rate  less  background  (made  on 
the  same  sample  using  the  same  problem,  system,  and  jeometry  numbers),  A.+.  is 
the  first  observation  made  after  observation  A.,  t  is  the  time  difference  between 
the  two  count  rate  observations,  and  X  is  the  effective  decay  constant.  It 
should  be  noted  that  the  value  of  X  for  the  second  to  the  last  card,  i  =  n-1, 
is  also  used  for  the  last  card,  i  ■  n.  The  time  difference  is  obtained  by 
multiplying  the  difference  in  days  by  24  and  adding  that  number  to  the  dif¬ 
ference  In  hours.  The  result  is  multiplied  b>  60  and  added  to  the  remaining 
difference  In  fractional  minutes.  The  effective  decay  constant  thus  determined 
is  used  only  for  computations  within  this  program  and  is  not  otherwise  trans¬ 
mitted  to  the  ofctput  data. 

Both  forms  of  DCISON  make  additional  tests  and  computations  to  correct 
for  decay  during  counting  and  to  determine  if  there  is  significant  decay 
between  observations.  It  is  arbitrarily  chosen  that  if  the  product  of  Xt  , 
tc  being  the  duration  of  count  observations,  is  greater  than  0.01,  the  ob¬ 
served  counts  (less  background)  are  corrected  for  decay  during  observation 
by  the  relation 


AXt 


(I 


-Xt* 

•e  ) 


where  A  is  the  observed  counts  and  A  is  the  corrected  counts.  For  each 
such  card,  the  observation  number  is  set  equal  to  one.  If  the  count  ob¬ 
servation  includes  background,  then  it  is  necessary  to  convert  sample  and 
background  counts  to  rates,  subtracting  the  latter  from  the  former,  before 
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making  the  above  correction.  If  the  product  of  the  time  difference,  t^, 
between  observations  and  the  decay  constant,  X,  Is  greater  than  0«.01, 
significant  decay  Is  said  to  have  occurred  between  observations.  If  such  I 
Is  the  case,  the  observation  number  Is  set  equal  to  one  and  the  time  at 
which  the  observation  Is  made  Is  unadjusted.  If,  however,  the  product 
Xtj  Is  less  than  0.01,  then  the  time  of  observation  Is  set  equal  to  that 
time  of  the  first  data  and  the  observation  number  Is  set  equal  to  that 
of  one  plus  the  data  Immediately  preceding.  The  source  programs  are 
listed  In  Table  XVI,  XVIII  and  the  output  listing  for  both  forms  of 
DC  ISON  Is  given  In  Table  XVI  lab,  XIX. 

Chi  Square 

CHISQ  Is  a  program  used  to  determine  If  a  counting  system  Is  performing 
satisfactorily.  This  program  obtains  Its  Input  data  from  DCISON  but  It 
Is  not  part  of  normal  data  reduction;  rather,  it  is  used  when  testing  a 
counting  system.  Data  for  CHISQ  is  obtained  as  a  series  of  count  rate  ob¬ 
servations,  each  for  an  equal  Interval  of  time.  It  should  be  noted  that 
any  counting  system  should  exhibit  only  the  Poisson  variations  of  random 
disintegration  process  except  for  minor  long  term  changes.  The  source 
program  is  listed  in  Table  XX,  and  in  Table  XXI  is  the  output  listing  for 
a  sample  program. 

Decay 

DECAY  is  a  program  used  to  correct  for  radioactive  decay  in  the  count l r. 
rate  of  a  sample.  In  the  present  scheme  of  analysis  three  programs  are  pro 
vlded  depending  upon  whether  there  are  one,  two  or  three  components  of  deca 
The  number  of  decay  components  ( i .  e . ,  decay  constants),  is  sensed  by  the 
preceding  program  DCISON,  and  appropriate  on-line  instructions  are  provided 
to  the  computer  operator.  For  this  program  the  cards  associated  with  the 
problem  constants  provide  the  decay  constants.  The  program  interprets  thes 
decay  constants  as  the  slope  of  straight  lines  on  simllog  paper  and  fits  th 
observed  counting  rate  data  to  the  straight  lines  such  that  the  sum  of  the 
straight  line  components,  one  for  each  decay  component.  Is  equal  to  the 
observed  counting  data;  as  a  function  of  time  the  smooth  curve  through  the 
counting  rate  data  is  optimized  by  the  least-squares  process.  The  source 
programs  are  listed  in  Tables  XXII,  XXIV,  and  XXXI,  for  Decay  I,  II,  and  il 
corresponding  to  I ,  2,  and  3  components  of  radioactivity.  The  program  es¬ 
timates  the  most  probable  counting  rate  for  each  component  of  the  radioacti 
sample,  up  to  a  maximum  of  three  components.  This  estimate  is  at  either 
a  given  day  and  time,  which  must  be  less  than  the  lowest  day  and  time,  or 
the  lowest  day  and  time  depending  upon  the  information  given  in  the  data 
of  priblem  constants.  Additionally,  each  data  point  is  corrected  back  to 
that  time  so  that  the  resulting  output  data  which  is  shown  in  Tables  XXIII, 
XXV,  XXVII,  respectively  for  Decay  I,  II,  d  III  consists  of  the  average 


counting  rate  for  each  component  as  well  as  the  average  total  counting  rate 
of  the  sample  along  with  the  counting  rate  of  each  observation  broken  Into 
Its  constituent  radioactive  components  along  with  the  sum  of  these  components. 

In  all  cases  the  other  experimental  parameters  such  as  time  ot  observation, 
background  counting,  and  time  of  observation  of  the  background  counts,  number 
of  samples,  and  number  of  observations  per  sample,  are  preserved  so  that  sub¬ 
sequent  programs  can  associate  confidence  limits  with  these  count-rate  es¬ 
timates.  Although  some  variance  of  data  points  from  the  Idealized  curve  Is 
possible,  extreme  variations  are  not  tolerated  by  the  present  program.  This 
Is  because  weighting  Is  not  employed  to  allow  counting  rates  that  more  closely 
approximate  the  predicted  value  (based  upon  known  half-life)  to  be  considered 
more  Important  than  '.hose  counting  rates  which  deviate  considerably  from  the 
predicted  values.  For  the  same  reason,  gross  errors  can  result  if  the  absolute 
magnitude  of  one  component  of  radioactive  decay  Is  small  In  comparison  with 
the  absolute  magnitude  of  the  sum  total  of  all  components  of  radioactive  decay 
present  In  the  sample.  These  limitations,  however,  do  not  impose  serious  res¬ 
trictions  for  practical  radiometric  analysis,  and  the  incorporation  of  these 
limitations  together  with  those  concerning  the  duration  of  time  of  observation 
and  number  of  observations  per  sample,  which  will  be  discussed  as  part  of  the 
following  program  VAR,  greatly  simplify  the  mathematical  analysis  and  re¬ 
quired  computer  programming. 

Analysis  of  Variance 

VAR  is  a  program  for  determining  the  average  counting  rate  of  a  series  of 
observations  made  over  a  period  of  time.  As  such,  there  are  possible  mechanical, 
sampling,  chemical,  (when  more  than  one  sample  is  being  compared)  and  non-Poisson 
counting  errors  in  addition  to  the  Poisson  variation  of  radioactive  decay  process. 
Its  input  is  the  count  rate  of  the  sample  and  background  counting  rate  pro¬ 
vided  by  DCISON.  The  program  also  estimates  the  confidence  interval  associated 
with  the  average  count  rate. 

The  program  considers  m  samples  whose  hypothetical  true  count  rate  is  the 
average  of  an  infinite  number  of  observations  made  on  an  Infinite  number  of 
samples.  Thus,  the  m  samples  represent  a  random  collection  from  the  sample 
population.  This  population  has  a  mean  and  a  standard  deviation,  whose 
estimates  based  ujgon  m  samples  observed  a  finite  number  of  times  are,  res¬ 
pectively,  p  and  a.  These  parameters,  along  with  identification  data,  are 
provided  as  output  data.  It  is  assumed  that  each  sample  is  observed  for  the 
same  number  of  times,  that  each  time  duration  of  observation  Is  the  same, 
and  that  a  separate  background  determination  is  made  for  each  observation. 

In  addition  to  determining  the  best  estimate  of  the  counting  rate,  an 
estimate  is  made  of  the  confidence  limit  associated  with  the  standard  de¬ 
viation,  0.90,  0.95,  and  0. 99-level  of  significance.  This  level  is  used 
in  calculating  the  upper  and  lower  limits  of  the  confidence  interval  associates 
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with  p.  Additional  ly,  the  number  of  additional  samples,  or  observations.  In 
the  case  of  single  sample,  required  to  reduce  the  estimated  error  of  the  measut 
ment  Is  calculated  and  punched  on  cards  together  with  the  other  output  data. 
The  program  Is  listed  in  Table  XXVI I  i,  for  the  case  of  one  sample,  and  the 
output  of  a  sample  program  Is  listed  in  Table  XXIX.  The  result  is  the  best 
estimate  of  the  counting  rate  of  each  component  of  radioactive  material,  up 
to  a  maximum  of  three,  at  a  given  day  and  time,  along  with  the  confidence 
interval  of  these  estimates.  This  best  estimate  Is  based  upon  a  series  of 
observations  made  on  a  series  of  samples.  Similar  source  programs  for  the 
case  of  many  samples,  are  listed  in  Tables  XXX  and  XXXII.  The  corresponding 
sample  output  data  listings  are  in  Tables  XXXI  and  XXXIII. 

Backing  Thickness 

BACKTH  calculates  the  thickness,  XF,  of  backing  material  for  radioactive 
samples  by  three  different  methods.  First,  the  data  may  already  be  given  in 
the  desired  units,  mg/cm  ,  which  merely  entails  punching  the  sample  number 
with  its  associated  thickness.  Second,  the  thickness  may  be  given  in  mils, 

In  which  case  the  density,  p,  of  the  backing  material  must  accompany  the  value 
of  thickness.  Third,  the  data  may  be  given  as  the  ratio  of  counts  observed 
when  the  backing  is  used  as  a  absorber.  Input  data  is  obtained  from  the  out¬ 
put  of  DECAY  I,  sinpe  this  program  is  also  not  part  of  the  regular  reduction 
routine.  If  the  density  is  missing,  tbis  indicates  that  the  backing  thickness 
Is  already  In  the  desired  units,  mg/cm  .  If  the  density  is  present,  then  the 
thickness  Is  In  mils,  which  requires  changing  these  units  to  milligrams  per 
square  centimeter.  This  is  accomplished  by  multiplying  the  given  value  of 
XF(l)  by 

10  ^ ( I n)  2.5b  (cm/ in)  p(mg/cm^) 

where  p  =  the  density  of  the  backing  material.  Thus,  XF ( I )  is  now  in  the 
desired  units.  On  the  other  hand,  if  the  value  of  XF ( I )  is  missing,  then  the 
data  card  must  contain  the  maximum  beta  energy  for  this  particular  isotope, 

E,  the  counts  observed  when  the  backing  is  used  as  an  absorber,  AA,  and  the 
counts  when  the  backing  is  not  used  as  an  absorber,  AO.  A  zero  value  of  XF ( I ) 
instructs  the  program  to  check  the  geometry  number  given  for  the  present 
execution  of  BACKTH.  A  geometry  code  number  less  than  10  indicates  narrow 
geometry.  If  this  is  the  case,  the  above  entries  on  the  current  data  card 
are  used  In  obtaining  a  value  of  the  mass  absorption  coefficient.  Of,  by 
employing  the  following  equation  for  narrow  geometry 

22  2 
cr  *»  >  cm  /gm 

where  E  =  energy,  Mev.  If  the  geometry  code  number  Is  greater  than  or  equal 
to  10,  indicating  wide  geometry  (2ir  or  bn) ,  then  the  given  value  of  the 
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maximum  beta  energy,  E,  ts  checked  to  see  ff  It  Is  less  than  1.45  Mev.  If  so, 
the  following  equation  for  wide  geometry  Is  used  to  determine  a 


a 


60 

e1.54 


If  the  value  of  E  Is  greater  or  equal  to  1.45  Mev.,  the  equation,  also  for 
wide  geometry,  Is  employed 


a 


200 

e6.25 


After  the  appropriate  value  of  a  ts  determined.  It  Is  used  In  calculating 
the  corresponding  value  of  XF(l).  This  is  done  by  substituting  the  values 
of  AO,  AA  and  a  Into  the  following  expression 


XF(|)  «  Jn-lAQi  _  .In  (AA) 


Immediately,  this  value  of  XF ( I )  and  Its  associated  sample  number  are  punched 
on  cards.  The  source  program  Is  listed  in  Table  XXXIV  and  Table  XXXV  lists 
the  output  data  for  a  sample  problem.  It  should  be  noted  that  these  programs 
do  not  handle  the  case  where  the  geometry  is  actually  neither  narrow  or  wide 
but  somewhere  between  the  two  special  cases. 

RESULTS 


The  response  of  the  selected  foils  to  moderated  nsutrons  is  shown  in 
Figure  4.  In  this  Figure  the  response  of  each  foil  was  normalized  by 
dividing  the  observed  activity  by  the  maximum  observed  activity.  This  normal¬ 
ized  activity  is  the  linear  left  ordinate  and  Is  plotted  versur  "he  paraffin 
thickness,  measured  in  inches  along  the  logarithmic  abscissa.  The  only  ex¬ 
ception  is  the  cadmium  ratio  determinations.  The  cadmium  ratio  w*s  calculated 
by  dividing  the  activity  of  the  bare  foil  by  the  activity  of  the  foil  covered 
with  20  mils  of  ctdmium.  This  ratio  is  the  linear  right  ordinate  plotted 
against  the  same  abscissa.  The  reciprocal  of  the  cadmium  ratio  Is  the  fraction 
of  the  total  number  of  neutrons  having  energies  greater  than  the  effective 
cadmium  cutoff  energy(N.B.,  for  the  cadmium  thickness  specified  above  and  for 
an  ambient  temperature  of  23°C,  the  effective  cadmium  cutoff  energy  is  0.48-  ev)  . 
Fporn  the  cadmium  ratio  an  estimate  may  be  visually  obtained  of  the  average 
neutron  energy.  This  may  be  used  as  a  comparison  basis  for  the  foil  response. 
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It  will  be  seen  In  Figure  3  that  the  response  of  each  foil  Is  fairly  narrow 
and  does  not  overlap  more  than  a  portion  of  the  adjacent  foil.  Further,  It 
will  be  seen  that  there  Is  no  energy  at  which  some  foil  does  not  have  a  sig¬ 
nificant  response.  Host  of  the  foils  have  responses  between  I  and  16  Mev. 

This  Is  the  energy  region  where  most  resolution  Is  needed  In  the  study  of 
the  penetration  of  W-Mev  neutrons  through  shielded  structures.-* 

CONCLUSIONS 

Twelve  feet  of  paraffin  will  thermal Ize  at  least  99  percent  of  all  14-Mev 
neutrons  as  evidenced  by  the  measured  cadmium  ratio-  of  Indium.  The  response 
of  each  selected  foil,  ordered  In  decreasing  energy,  overlap  no  more  than  the 
response  of  the  adjacent  foils  and  there  Is  no  energy  "gap"  for  which  there 
Is  no  response  of  at  least  one  foil.  The  resolution  of  foil  response,  I.e., 
the  width  of  the  response  curve.  Is  most  suitable  for  measuring  neutron 
streaming  through  ducted  entranceways  of  protective  structures. 

A  series  of  computer  programs  Is  most  useful  for  the  reduction  of  foil 
activation  data.  The  programs  are  written  for  particular  application  to 
shielding  studies  where  several  different  foils  are  irradiated  in  many  dif¬ 
ferent  locations  In  the  shield.  The  programs  take  care  of  the  "bookkeeping" 
for  such  an  experimental  program,  such  as  foil  positions,  foil  type,  radio- 
metric  system  and  geometry,  and  different  observations  for  the  statistical 
analysis  Incorporated  In  the  program. 
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ORGANISES  DATA  BY  PROBLEM, 
CHANGES  OATA  format 


[_RAW 

DATA 

PRE  DATA  COMP 

POST  DATA  COMP 

DATORG  1 

DATORG  II 

DATORG  III 

RESOLUTION 

OROERS  BT  SAMPLE 


ORDERS  BT  SYSTEM 


ASSOCIATES  PROBLEM  INFORMATION 
WITH  PROBLEM 


ASSOCIATES  SYSTEM  INFORMATION 
WITH  SYSTEM. 


CORRECTS  FOR  TOTAL  COUNT  LOSS 


NO 


DCISON  II 

SUPPLIES  BAOHOROUN 
COlWTS  where  MISSI 
DECIDES  IF  DECAY  iS 
SIGNIFICANT  ESTlMATi 
A  OCCAY  CONSTANT 


DECAY  II 


DECAY  III 


PREP 

VAR  II 

VAR 

II 

gets  'jslues  from 
tab,  sets  up  data 

FOR  VAR  II 

computes  a  mean 
counting  batf  a 
CONFIDENCE  LIMITS, 
OECIDFS  IF  DATA  IS 
POISSONLY  OiSTRIBUT 


DETERMINES 

thickness  of 

BACKING 

material 


COMPUTES  A  MEAN 

count  ng  rate,  a 

CONFIDENCE  LIMIT 


Figure  i.  Pronr.im  Serjin  ncr  . 
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gure  Z.  Format  of  Counting  Input  Data. 


COUNTING  DATA 


NOTE 


NOTE  AS  MANY  CAROS  AS  T HERE 
ARC  2[S  IN  COLUMN  2  IN 
PRCCECOiNG  DECK 


AS  MANt  CANOS  AS 
PROBLEMS 


XXX 

xxxx 
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H-13 

LOW 

LOW 

PR  OB 

DAY 

TIME 

NO. 

NOTE 
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r€CK - — 


X 

COL  2 

XXX 
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»:  SEARCH  FOR  LOW 

. . 1 

TIME  , 

OECAr 

PR08 

SECONO  OAT 

cooc 

NO. 

2  *lOR  TIME  ANO  DAY 

GIVEN 

NOTE  AS  MANY  CANOS  AS  THE  TOTAL 
NUMBER  OF  GROUPS  EON  ALL 
PROBLEMS 


XX 

COL  4-5 
10  NO  OF 

EA  SAMPLE 


XXX 

SC 

RROB 

NO. 


NOTE  ONE  CANO 


NOTE  AS  MANY  CANOS  AS  THE  SUM 
OF  COL  5  IN  PRECEDING  SET 


XX 

COL  4-5 

NO  OF  SAMP 

IN  EA  GROUP 


XXX 
•  -IO 

PROS 

NO 


XX  XX 
CCL  2-5 

TOTAL  NO  OF  GROlPS  FOR 
ALL  PROBLEMS _ 


XXX 
:ol  2-4 

LOR 

DAY 


XXXX 

XXX 

XXXX 

XXXX, 

XX. XX 

XX 

1  B-4 

11-13 

»-!• 

21-25 

28-32 

34-40 

LOW 

HIGH 

HIGH 

BAG 

BAG 

SYSTEM 

TIME 

DAY 

TIME 

COUNTS 

TIME 

NO 

X 

|COL  5 

BACKGROUND 
CODE 


0-  DATE  GIVEN  WITH  COUNTING  DATA- 
X  x  lOCX  AT  FOLLOWING  CANOS,  VALUE 
OF  X  INCXCA7ES  NO  OF  BAG  CAROS 


XXX  X  .XXXEtXX  #XXXE!V 
COL  2-4  6  B-i5  17-24 

PROS  NO  OF  ' - DECAY  CONSTANTS 

NO  CCCAY  CONST  THIS  PROBLEM -a  MAXIMUM 


XXX 
3-5 
NO  OF 
PROBLEMS 


X.XXXXEtXX  XX 
COL  2-H 

Nf  SOLUTION  SYSTEM 

TIME  NO 


COL 


XX 
4-5 
NO  OF 
SYSTEMS 


NOTE  AS  MANY  CAROS 


CAROS  AS  SYSTE 


PROBLEMS 


NOTE  ONE  CARO 


NOTE  AS  MANY  CAROS  AS  THERE  ARC 
SYSTEMS. 


NOTE  ONr  caro 


Figure  3.  Format  of  problem  definition  cards. 


Paraffin  Thickness  (inches) 


Normalized  Foil  Response 


Figure  b.  Response  of  selected  foils  to  moderated  neutrons. 


TABLE  I.  List  of  Foils 


Foil  Reaction 

Dauqhter  Half-Life 

lnM5  (n,n‘)  In11* 

4.4  hr 

.n”5  (n,Y)  .nN6 

54  min 

Cu63  (n,2n)  Cu62 

9-9  mtn 

.,27  (n,2n)  ,n’26 

13.3  d 

S32  (n,p)  P32 

14.3  d 

S34  (n,a)  SI3’ 

2.62  hr 

P31  (n,p)  SI31 

2.62  hr 

Fe38  (n,p)  Mn38 

2.56  hr 

Mo32  (n,2n)  Mo^1 

15*5  min 

Al27  (n,p)  Mg27 

9.51  min 

Al27  (n,a)  Na24 

15  hr 

Al27  (n,y)  Al28 

2.3  min 

Cu65  (n  y)  Cu66 

5*1  min 

TABLE  I!  PRE  DATA  COMPILER  SOURCE  PROGRAM 


C  PRE-DATA  COMPILER 

DIMENSION  MM (125) 

X«l 

10  READ  1Q5  «N1 

105  FORMAT (15) 

PUNCH  105, N1  j 

DO  200  I  *  1 ,N1 

READ  106, NYX  , 

200  PUNCH  106, NYX  I 

106  FORMAT <40H 

132H  ,  I 

READ  105, N 
PUNCH  105, N 
DO  201  I a 1 ,N 
READ  807  »MM ( I J 

201  PUNCH  807, MM< I ) 

N3*0 

DO  202  1=1, N1 
READ  105, N2 
N3«N3+N2 

202  PUNCH  105, N2 
DO  203  I  *  1 , N3 
READ  1 06 , N YX 

203  PUNCH  106, NYX 

READ  105,  N5  I 

PUNCH  105  *N5 

N6*0 

DO  205  1*1, N5 
READ  40 1 *N7 »NNM 
N6*N6+N7 

205  PUNCH  40 1 ,N7 ,NNM 
DO  206  I s 1 ,N6 
READ  106, NYX 

206  PUNCH  106, NYX 

JmQ 

DO  208  I  *  1  »N 

READ  408 , IC  ,  I P 

IF ( IC-2)208,209,208 

209  U«J+1 

408  FORMAT (12,14) 

208  PUNCH  4Q8»IC,IP 
DO  210  I»l#J 
READ  106, NYX 

210  PUNCH  106, NYX 
IZ«-1 


PUNCH  109.12 

109  FORMAT (40X.32X.I2 ) 

3  READ  100.  B4.B3.1A7.IA6.IA2.IA5.D.B6.IB5.IA3.IA4.1A1.IB2.B1 
IF!  IB5)5.30.5 

30  PUNCH  110 

110  FORMAT I40X.40X) 

PRINT  808 

808  FORMAT (27HEND  OF  PROBLEM.  REREAD  DATA  J 
PAUSE 
R-IC+1 

IF (fc-N) 10.10.300 
5  IFUB5-MM<Kn3.4.3 

4  C2-IB2 
B2-C2/100. 

IF ! B4 )  1.2.1 

1  PUNCH  101.  IA1.IA2.IA3. I A4. IA5.IA6.IA7.B1.B2.B3.B4 .185. D.B6 
GO  TO  3 

2  PUNCH  102,  IA1.IA2,IA3,IA4,IA5,IA6,IA7,B1,B2.IB5,0,B6 
GO  TO  3 

300  IZ--9 

PUNCH  109.12 
PRINT  500 
STOP 

500  FORMAT ! 13HLOAD  POST  PDC ) 

100  FORMAT(F7.2,F8.0,3I3, 1 2 ,E9.3 ,F5 .2 ,2 1 4 , 1 5 , 14 , 14, IX ,F8 . 0 ) 

101  FORMAT! 14, 1 3 , 1 4 , I  5 , I  2 ,2 1 3 ,F8 .0 » F7 .2 , F8 . 0 . F7 .2 . 1 4 , E9 . 3 , F5 .2 > 

102  FORMAT! 14 . 1 3 . 1 4 , I  5 . I  2 .2  I  3 ,F8.0 .F7. 2 . 1 5X  ,  1 4 ♦ £9 .3 ,F5. 2) 

401  FORMAT ! 2 1  5 ) 

807  FORMA T ! I  4. 35H  ) 

END 
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TABLE.  IV  POST  DATA  COMPILER  SOURCE  PROGRAM 


C  POST  DATA  COMPILER 

DIMENSION  NSUOO> 

2  1*2 

l  READ  lOOtNYX 
PUNCH  100 .NYX 
IF  INYX+l  U.10.1 

20  READ  101*  IAl .IA2.fA3.IA4.lA5.IA6.IA7.Bl *B2.B3,B4*IB5*D*B6.NYX 
NS  I  H-IAl 

11  READ  101*  tAl.IA2.lA3.IA4,IA5»IA6,!A7»81,B?,B3.B4,JB5.D.B6.NrX 
IF i I A 1 i 16*15*16 

16  U*i 

13  IFlNS' JJ-TA1U2.11.12 

12  J* J*1 

I F  *  J- I » 13.13.14 

14  !  =  I ♦  1 

NS  J I >  =  ! A 1 
GO  TO  11 

15  PRINT  102 

I*  (  1-1  ><*00. 400. 408 
408  N*l 

DO  300  I  =  2  *  N 
J»l 

402  IFlNSt J)-NS(J-l) 1401*300.300 
401  AEP->NS<  j  » 

NS  <  J 1 «NS <  J- 1 ) 

NSIJ-l i»XEP 
IF ( J-2 140.300.40 
40  J-J-l 

GO  TO  *02 
300  CONTINUE 
1-N 

400  DO  2000  J*1*I 

200  READ  100. NYX 

If  <NYX. 1 ) 200.204,200 

20'-*  READ  101*  IAl.IA2.IA3.1A4,lA5.IA6*lA7,Bl,B2*B3*B4,IB5.D»B6,Nvx 
IF  !  IAU202.2001.202 
202  1*  !  IA1-NSI J)>204« 203.204 

20 1  PUNCH  10’ ,1A! ,IA2*IA3*IA4,IA5,1A6.IA7.B1 ,B2,B3*B4,TB5.D.B6.NYX 
GO  TO  204 

2001  IF  iO-I ) 2002 .2000. 2002 

2002  PRINT  102 
■  iOOO  CONTINUE 

PRINT  103 

Punch  101*IAI*IA2»IA3*IA4»IA5*IA6.IA?»81«B2*63*B4*I85*D.86»NYX 
PAUSE 


GO  TO  2 

100  FORMAT »40H  , 

132H  *12} 

101  FORMAT ( 14,I3,|4,I5.I2,2I3.F8*0.F7.2.F8.0.F7.2*I4,£9.3,F5.2. 

102  FORMAT ( 1 1HREREAD  DATA) 

103  FORMAT < 16HREAD  NEW  PROBLEM) 

END 
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C  DATA  ORG  II 

1  READ  200  »NYX 
IF(NYX+1 1 1  •  2 • 1 

2  READ  201 • IBS 
IX5*IB5 
PRINT  202 
PAUSE 

READ  100  »N1 
PUNCH  100. N1 
DO  10  I»1*N1 
READ  200. NYX 
10  PUNCH  200.NYX 
READ  100. N2 
DO  20  1=1 *N2 

READ  20<t  «  I °ROB .NLM.XLAMl.XL  AM?  «  XL  AM? 
I F  ( I PROB- 185)20*21.20 

21  IY1=IPR0B 
I Y2=NLM 
XX1=XLAM1 
XX2=XLAM2 
XX3«XLAM3 

20  CONTINUE 
NN*0 

DO  22  1*1. N1 
READ  2Q5.N3.N4 
PUNCH  205.N3.N4 

22  NN*NN+N3 

DO  23  1=1* NN 
READ  200.NYX 

23  PUNCH  200 »NYX 
NYX--2 

PUNCH  203  »NYX 

PUNCH  204. I Y1 *IY2.XX1  .XX2.XX3 
READ  100.N5 
NM*  0 

DO  25  1  =  1. NS 
READ  2Q5.N3.N“ 

NN  *NN>N2 

IF  <N4-Ia5»25. 26.25 
2b  PUNCH  205»N3»N4 
25  CONTINUE 

DO  27  1=1. NN 
READ  205.N3.N4 
IF (N4-IX5  »27.?8.2  7 
28  PUNCH  205.N3.N4 


27  CONTINUE 
NYX»**3 

PUNCH  203  »NYX 
NN-0 

DO  29  1-1.N2 
READ  209  ,N 
IF (N-2)29»30»29 

30  NN*NN+ I 
29  CONTINUE 

»F <NN)31 *32*31 

31  DO  33  I«1 »NN 
READ  206  »  J1 ,  J2  * J3 
Ip l  .,3-1X5 J33.34.33 

34  PUNCH  206.J1.J2.J3 
33  CONTINUE 

32  READ  200. NYX 
PUNCH  200. NYX 
IF(NYX+1)32. 35.32 

35  READ  207.IA1.IA2.IA3.IA4, IA5.IA6.IA7.B1. B2.B3.B4.IB5.DK. B6 
I F (  IA1J36.37.36 

36  PUNCH  207,  I A1 • I A 2 , I  A  3 , I A 4 , I A  5 ♦ I A6 , I A 7 , B 1 . B2 .83 ♦ 64 . 1 8 5  ,DK .B6 
GO  TO  35 

37  NYX«*-9 
PUNCH  203. NYX 
PRINT  203 
PAUSE 

GO  TO  1 

100  FORMA  T ( I5> 

200  FORMAT  { 40H  . 

1 32H  ,12) 

201  FORMAT (40X.14X.I4) 

202  FORMAT ( 11HREREAD  DATA) 

203  FORMA  Tt40X.32X.I2) 

204  f ORMAT ( 14. 12  .3E9.3 ) 

205  F  ORMA  T ! 2  I  5 ) 

206  FORMAT (  14. 15. 14) 

20  7  FORMAT  t  l^,I3.I4.I5,I2,?I3,F8.0.F7.2,F8.O,F7.2,I‘*.E9.3.F5.Z) 

208  FORMA  T  U4HL.CAD  NEXT  PR 08) 

209  FORMA  I  I  J  2  ) 

END 
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TABLE  X  DATORG  III  SOURCE  PROGRAM 


C  DATA  ORG  III 

DIMENSION  TAUI 50 ) ♦  I  SYS  I  50 ) .INTO! 100) , I  NTT ( 1 00 ) ♦ I  PD  11 00 > , BKGC f 1 00 ) 
DIMENSION  BKGT ( 100) * IFT(IOO) *  I  SMI  100  I 

1  READ  100 .N 1 *NYX 

2  DO  1C  1*1. N1 

10  READ  101 tTAU(I) .ISVSI I) 

NX-0 
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11  NX*NX+N2 
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15  READ  103.NYX 
PUNCH  103. NYX 

I  F *  NYX  + 1 >15.14.15 
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34  I X2= I A2 
I  Y5= 1 85 
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16  CONTINUE 
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50  CONTINUE 
NYX»-5 
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53  IFINYX+9J60, 54.60 
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PUNCH  300. NYX 
GO  TO  34 

54  PUNCH  300. NYX 

300  FORMA T (40X.32X.I2 ) 
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204  FORMAT { 21HEND  DATA  ORG  PROGRAMS ) 
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TABLE  XII  RESOLUTION  SOURCE  PROGRAM 


C  RESOLUTION 

1  READ  201 »NYX 
PUNCH  201 *NYX 
IF<NYX+1)1»10.1 

10  READ  lOOtTAU.lSY 

11  READ  201 #NYX 
PUNCH  20 1 »NYX 

IF  (NYX+5 )  11 *12  *11 

12  READ  101.IA1,IA2.IA3.IA4,IA5«IA6,IA7,  Bl*B2*B3*B4»lB5«D«B6,NYX 
IF « 1SY-IA2I13. 14*13 

13  PRINT  104 
PAUSE 

GO  TO  1 

14  C*B1/B2 
FTAU«1«“C*TAU 
TRUE«=Bl/FTAU 

PUNCH  lOltlAl , IA2  * IA3  ,  IA4,  IA5, IA6,  IA7,TRU£%B2,B3.B4,IB5,D,B6,NYX 
READ  l01tIAl*IA2*IA3»IA4,IA5*IA6«!A7,  BltB2«B3»B4*lB5*D»86»NYX 
1F<NYX+6)16. 15*16 
16  IF ( NYX+9 ) 14»20»14 
20  PRINT  105 

PUNCH  202.NYX 

PAUSE 

GO  TO  1 

15  PUNCH  20P.NYX 
GO  TO  10 

100  FORMAT (Ell«4»I4l 

101  FORMAT ( I4,I3»I4,I5, I  2 . 2  I  3 , F8 • 0 »F 7 . 2 , F8 .0 .F7 . 2 » 1 4 . E9* 3 . F5 . 2  *  I  2 ) 

104  FORMAT ( 27HERROR.  MIXED  SYSTEM  NUMBERS) 

105  FORMAT « 16HL0AD  NEW  PROBLEM) 

201  FORMAT ( 40H  » 

132H  ,12) 

202  FORMAT  <40X,32Xtl2l 
END 


r  ABLE  XIII  OUTPUT  RESOLUTION  -  INPUT  DECISION 
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TABLE  XIV  DEAD  SOURCE  PROGRAM 


C  GM  DEAD  T I ME  LOSS  CALCULATION 

DIMENSION  AMP (200) .AN l 200 ) ♦ I T ' 200 > 

90  READ  100. NYX 

READ  107.IPRB.NRLAM.XLAM 
107  FORMAT (I4.12.E9.3) 

10  READ  100. NYX 

I F  <  NYX  1 1 ) 10.11.10 

11  READ  101.1A1.IA2.IA3.IA4.IA5.IA6.IA7.B1.B2.83.BL.NYX 
I F  (  I  A5-*7 }  71 .12  .71 

12  1 SMP- I A 1 
l SYS* I A2 
1  DY^»  i  A  3 
IGEO* ! A6 

I  »0 

GO  TO  200 

201  READ  101.IA1.IA2.IA3.IA4.IA5.IA6.IA7.B1  »B2.B3.Bl»NYjc 
IF ' I A 1  —  I SMP  *205*202.205 

202  IF< (Svs-IA2I205.203.20S 
2C3  If  UDY-IA3)  205.20L.205 
?CA  IF( IGEO- l A6 ) 205 .200 , 205 
200  IU*1 

AMP  I ! ) »B 1 
222  LFT*lA<r/10Q 

ITU  |  «IA4-LFT*100+LF  T+60 
GO  TO  201 
71  PRINT  103 

103  FORMAT i 12HNOT  TAU  DATA) 

PAUSE 
GO  TO  90 
205  J»1 

AN ( J ) *AMP i J ) 

N*J-1 

DO  200  I-l.N  < 

TIME- IT( J)-IT(  I  ) 

300  ANU ) -AN) U»*EXP(XLAM*TIME) 

SUMT  «*  0.0 
SUMB  «  0.0 
DO  *0  I-l.J 

ST  a  AMP ( l ) * ( 1 *UAMP U ) /AN ( I ) ) 

SUMT  •  SUMT  t  SI  ♦ 

S6  »•  AMP  (  l  1**2 
*♦0  SUMB  *  SUMB  +  SB 
DT  *  SUMT  /.SUMB 
IF  <DT  >400.402*402 
LOO  DT  =0 * 


402  PUNCH  106 

PUNCH  105.ISYS.IDY.DT 
13  IF(NYX+9) 700,70,700 
70  PRINT  102 

102  FORMAT ( 16HREAD  NEW  PROBLEM) 

60  TO  90 

700  IF(NYX+6)7Q1»11»7Q1 

701  PUNCH  100. NYX 
60  TO  11 

100  FORMAT {40H  • 

1 32H  .12) 

101  FORMAT < I4,I3,I4,I5,!2,2I3.F8.0,F7.2»F8.0,F7.2,18X.I2) 

105  FORMAT (2I5.2X.E10.4) 

106  FORMAT ( 6HSYSTEM.4H  DAY.9H 
END 


TAU) 


TABLE  XV  DEAD  OUT 
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*8LE  XV!  DECISION  I  SOURCE  PROGRAM 


DECISION  I 
DIMENSION  X(  30.6) 

202  FORMAT ( 16HLOAD  DECISION  II) 

203  FORMAT ( 16HBKG  INFO  MISSING) 

112  FORMAT ( 16HREAD  NEW  PR08LEM ) 

300  READ  200. NYX 

READ  201.PN.NRLM,  XL.XLAM2.XLAM3 
IF  (NRLM-1  )250»350,250 
250  PRINT  202 
STOP 

350  PUNCH  200.NYX 

PUNCH  201.PN.NRLM.  XL  .  XL AM 2 . XL AM 3 

210  READ  200. NYX 
PUNCH  200. NYX 

IF  INYX  +  1 ) 210.211,210 

211  L  L  *0 
KKK«0 

221  READ  100,ICl,IC2,IC3,IC4.C5,C6,NYX 
IF(NYX+5)212. 3^,212 

212  I =LL ♦  1 
X(I.1)»IC1 
X(  I,2)*1C2 
X( 1 , 3 | « I C  3 
XII  .M»!C4 
X ( I .5 ) *C5 
X(I,6)»C6 
LL*LL+ 1 

GO  TO  221 

34  READ  102,  IA1,IA2,IA3,IA4,IA5,IA6,IA7,B1,B2,B3.B4,IB5,XLL,B6,NYX 
IF (XKK) 4000,4004,4000 
4004  KICK  =  1 
4001  I G* I A6 
LRRa  I A5 
I SM= I  A 1 
ID“0 
I  T«0 
Jol 

4000  IF (NYX  +  9)215, 216,215 
216  PRINT  112 

PUNCH  877, NYX 
877  FORMAT ( 40X  ,  32X  ,  I  2  ) 

IF(LP.R-7)7001, 7002.7001 
7002  PPINT  7004 

7004  FORMAT! 17HTHIS  PROB  TO  DEAD;  , 

7001  *>AUSE 


CO  TC  300 

215  IF  'NYX*6>4002*270,4002 
<♦002  IF  UA6- IG1  4001  »fc  »<*001 

6  IF UA1-ISM)4001*53»4001 
270  PUNCH  8  7  7  * N  V X 

GO  TO  211 
53  IF  I  ID)  7  *46,  7 
46  101= IA3 
i r l- i aa 

7  IFUA5-1)  8,10,13 

13  IFUA5-2J  8,10,9 

8  PRINT  104 
PAUSE 

GO  TO  34 

9  I F  < B3 )  11,12.11 

11  IF(B4*  22,12,22 

12  A3=*  I A  3 
A4sjA4 
1*1 

GO  TO  14 

17  I F ( l -LL  >  1 5  » 1 5  *  1 6 

15  IM*1 

14  IFIA3-XI 1,1)1  17,18,19 

18  I F ( A4-X  f  I » 2  )  )  17,19,19 

19  IF ( A3-X( ! ,3) )  20,20,17 

20  Ir(A4-XII, 4))  21,21,17 

16  PRINT  109 
GO  TO  10 

21  B3=X( 1,5) 

B4*X( 1,6) 

I A  3  *  A  3 
IA4-A4 

22  C=B1/B2-B3/B4 
GO  TO  64 

10  C-B1/B2 
64  U=XL*B2 

IHIJ-.Oli  28,29,29 

29  C'C*U/< 1 ,0-EXPl-U) ) 
IA7-1 

28  1 F ( I A3- ID  I  30,31,32 

30  PRINT  107 
GO  TO  32 

31  IF(JAfc-IT)  30,32,32 

32  IF  I  ID)  33,36,33 
36  A«=24»  (  I  A3*- 1 D  ) 

B-»  I A4-  l  T 
GO  TO  37 

33  Ar24#< IA3-ID1 ) 


B-IA4-IT1 

37  D»0. 

K-l 

GO  TO  38 

41  D«D+100. 

K*K+1 

38  IF(B-D)  39*40.41 
40  BO-O. 

GO  TO  42 
61  D*D-10G. 

67  BD*B-D 
GO  TO  42 

39  IF(K-l)  60.60.61 
60  D=-100. 

65  IFfB-O)  66,40.67 

66  0*0-100. 

GO  TO  65 

42  T=60,*A+.6*D+BD 
IF<XL*T-.01)  43,44,44 

44  1A7*1 

GO  TO  47 

43  I  A7  =  J 
J*J+ 1 

49  PUNCH  102»IA1»IA2,ID1»IT1,IA5,IA6,IA7,  C.82 ,B3 ,B4 , 105 .XL  ,B6 .NYX 
ID*IA3 
IT* I A4 
GO  TO  34 

47  PUNCH  102.IA1,IA2.IA3,IA4,IA5,IA6,1A7,  C ,B2 ,B3 ,B4 , IB5 , XL  ,B6 .NYX 
ID* IA3 
I T* I A4 
101= IA3 
I T 1= I A4 
GO  TO  34 

100  FORMAT ( 14,I5,I4,I5,F8.0*F7.2,39X,I2) 

102  FORMAT ( 14,13.14,15, 12,2 I3.F8.0.F7.2.F8.0.F7. 2. I4.E9.3.F5. 2. 12) 
104  FORMAT  <  28HERROR ,  BKG  SUBT  DATA  MISSING) 

107  FORMAT (30HERROR,  INCORRECT  TIKE  SEQUENCE) 

109  FORMAT (46HERR0R,  DATA  CARD  DOES  NOT  FIT  GIVEN  TIME  RANGE) 

201  FORMAT ( I4.I2.3E9.3) 

200  FORMAT ( 40H  , 

1 32H  ,12) 

END 
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TABLE  XVI H  DECISION  II  SOURCE  PROGRAM 


C  DECISION  II 

DIMENSION  XI  30.6) 

200  FORMAT <40H 

132H  .12) 

20 2  FORMAT < 15HL0A0  DECISION  I) 

203  FORMAT ( 16H8KG  INFO  MISSING) 

112  FORMAT (16HREAD  NEW  PROBLEM) 

NY-0 

300  READ  200.NYX 
PUNCH  200. NYX 

READ  201 .PN.NRLM.XLAM1 * XLAM2  . XLAM3 
PUNCH  201 *PN  »NRLM  »XL AMI .XL AM 2. XL AM 3 
I F ( NRLM- 1)210.250.210 
250  PRINT  202 
STOP 

210  READ  200. NYX 
PUNCH  200 .NYX 

IF INYX+1 )210.211.210 

211  LL-0 

I  Yl  =  0 

ID-0 

IT-0 

K-2 

J*1 

221  READ  100.IC1.IC2.IC3.IC4.C5.C6 .NYX 
IF<NYX+5)212»34,212 

212  I-LL+1 
X(I.1)-IC1 
XI I .2 )-IC2 
X ( I . 3 ) • IC3 
XI I »4)«IC4 
X ( I .5 ) -C5 
XI  I . 6 ) -C6 
LL-LL^l 

GO  TO  221 
83  I Yl= I A1 
IY2-IA2 
I Y3= I A3 
I Y4= IA4 
I Y5= I A5 
IY6±IA6 
IY7-IA7 
Gl-C 
Y2-B2 
Y3-B3 


Y4»B4 
JY5* IBS 
Y6=B6 

34  READ  102*  I Al . I A2 * IA3 • 1 A4. I A5 • I A6. IA7.81 .62 *B3 *84 , 1B3 . 0.B6 .NYX 
I F  C IA1)53.1000.53 
1000  K  =  3 

GO  TO  84 

701  PUNCH  U5.NYX 
115  FORMAT (40X.32X.I2) 

IF (NYX +9)211. 216.211 
216  PRINT  112 
PAUSE 
GO  TO  300 
53  IF(!D)7,46«7 
46  ID1*=IA3 
I  T 1  =  I A4 

7  IF  (  I  A5**l  )  8.10.13 

13  I F ( I A5-2  1  8.10.9 

8  PRINT  104 
PAUSE 

GO  TO  34 

9  I F ( B3 )  11.12.11 

11  IF<B4)  22.12.22 

12  A3MA3 
A4*  I A4 
1*1 

GO  TO  14 

17  IF(  I-LU15.15.16 

15  I ® I 4 1 

14  IF ( A3“X I  I .1) I  17.18.19 

18  IFCA4-XI I  .2) »  17.19.19 

19  IF ( A3-X I  1.3) )  20.20,17 

20  I F ( A4“X ( I .4) |  21.21.17 

16  PRINT  109 
GO  TO  10 

21  B3*X (  I  *5) 

B4-X ( 1.6) 

IA3-A3 

I A4s A4 

22  C-81/B2-B3/B4 
GO  TO  23 

10  C=Bl/82 

23  IF  C K—2 I  1001.57.1001 
57  ft*l 

NMU*  0 
GO  TO  83 

1001  IF UY1-IA1 11003.1002. 1003 

1002  IF f  I  Y6-IA6)  1003.58.1003 


1003  K«2 

00  TO  84 

58  W«24*( IA3-IY3I 
Y-IA4-IY4 
Z-Q. 

j-i 

00  TO  59 

62  Z-Z+100. 

J*J+1 

59  IF ( Y-Z )  60*61*62 
61  YZ-O. 

00  TO  63 

91  Z«Z-100. 

94  YZ* Y“Z 

00  TO  63 

60  I F  C  J- 1 1  90*90*91 
90  Z*“100. 

92  IF  <  Y**Z  )  93.61*94 

93  Z-Z-100. 

00  TO  92 

63  T-60.*W+.6»Z+YZ 

XXL* LOG CGI ) /T-LOGCC ) /T 
84  V-XXL*Y2 

IFCV~0.1I  65*66*66 

66  G1*G1*V/C1.0~EXPC~V) ) 
J«1 

65  1FCIY3-10I  67*68*69 

67  PRINT  107 
PAUSE 

GO  TO  34 

68  IFCIY4-ITI  67*69.69 

69  IF ( ID)  70.71.70 

71  A*24*C IY3-ID) 

B-IY4MT 

GO  TO  72 

70  A»24*C IY3~ID1| 
B-IY4-IT1 

72  D«lC0. 

00  TO  73 

76  0>D+10G* 

73  IFCB-D)  74*75*76 
75  BD*Q h 

00  TO  77 

74  0*04100* 

B6-B-D 

77  T«6Q«*A*.6»D+9D 

IF  CXXL*T-.C1 )78»79»79 
79  IY7«1 


J 


83 


GO  TO  80 
’a  IV^sNHUrl 
NMU»NMU*1 

PUNCH  102*  ]rl,JY2.!Dl.ITl*iY5,IY6»IY7,Gl*Y2*Y3*Y4*jv54Xxt.  *Yc«NY 
ID*IY3 

II “IY4 
GO  TO  1009 

80  PUNCH  102*  IY1*IY2*IY3»IY4*IyS*!Y6*IY7*G1»y2*Y3*Y4*JYS,XXI..*Y6*NY 
NMU<=  1 
I0C I Y3 
IT--IY4 
I Dl= I Y  3 
JTl=lY4 

1009  IFHC-Ii  100"’*83*1007 

1007  IF(K-2)701*1 00 8*701 

1008  ID=0 
I  T  *=0 

GO  TO  53 

100  FORMAT ( 14* 15*14,1 5*F8.0*F 7*2* 39X* 12 ) 

102  FORMAT  (  1 4  *  1 3  *  14 *  I  5  *  I  2 *2  I  3 ,F8.0 *F7 . 2 ,F8 *0*F7 *2 » J 4  , E9* 3 .F 5 . 2 , I  2 > 

103  FORMAT { 17HERROR,  MIXED  DATA) 

104  FORMAT ( 28HERROR,  BKG  SUBT  DATA  MISSING} 

105  FORMAT ( I  4  »E9*3 i 

107  rORMA  T  <  30HERROR  «  INCORRECT  TIME  SEQUENCE} 

109  FORMAT (46HERR0R*  DATA  CARD  DOES  NOT  FIT  GWEN  TIME  RANGE) 

114  FORMAT (4SHUSE  SAME  DECAY  CONST  FOR  LAST  AND  2ND  TO  LAST , 1X4HCARD » 
201  FORMAT  < I4*I2*3E9*3) 

END 


TABLE  XIX  OUTPUT  DECISION  II  -  INPUT  DECAY  II-III 
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table:  XX  CHISO  SOURCE  program 


C  CH!  SQUARE  SECTION 

DIMENSION  TOO;  »SN « 30 ) ,B < 4 » ,C ( 4 » 

B;l J-21.S 
B 1 2 ( ! 14*6 
B 1 3 1  * 4,00 
B  i  4  k  =  2,00 
C  5 1  »  sao.O 
C  *  2  »  *27*2 
C  k  3 ) - 1 1 «6 
C  (4 » -=7*50 
1*1 

302  READ  205  ,NY* 

KK  *  1 

READ  201 ,PN,NRLM,XLAM1 ,XLAM2 ,XlAM3 
*  IF (NRLM- 1 »301 ,300,301 
301  PR3  NT  304 

304  EORMAT;30HTHIS  DATA  WITH  DECAY  II  OR  lift 
PAUSE 

GO  TO  302 
300  READ  205  »NYX 

IFINYX+1 ,300,305,300 

305  READ  800,  JAl«lA2,IA3* 1A4,IA5,IA7,IA8*B1 .B2 ,83 »B4 , \ B5 ,D *86 ,N YX 
GO  TGI  102*100*102  » 100 ) * IA5 

100  PRINT  309 

309  FORMAT  I 27HNOT  CHISO  DATA, LOAD  DECAY  I) 

PAUSE 
GO  TO  302 

80  READ  800,  IA1»IA2,  IA3»  IA4,IA5,1A7,  IA8,B1  ,B?,B3,B4, 1 95  ,D ,3*>  ,N  YY 
IF (NYX*9i311, 315,311 

311  IF (NYX+6»313«316*313 
313  IFMC.EP-*IA1)12,81,12 

81  IF.KEP1-IA2)  12,82,12 

82  I F  ( XEP2- I  A3 )  12,83,12 

83  IF<KcP3-lA4>  12,84,12 

84  IFIKEP4-IA7I  12,85,12 

85  I F  4  KEP 5 - 1 85  »  12,1,12 
1  $Nm=Bl 

T « I  M-B2 
1  «1  +  1 
GO  TO  80 
102  KEP*1A1 
REP1 *  1 A2 
KEP2* I A  3 
KEP3=IA4 
*CEP4=  I A7 


fcEP5=IBS 
GO  TO  1 

12  MA= 1-1 

IF (MA-2J500.13.13 

13  DO  40  I *2  »MA 

IF  ( T  ( 1  J-TlIJ  J41.40.41 
41  PRINT  104 
GO  TO  102 
40  CONTINUE 
500  MT  =MA/ 10  +  1 
DO  70  I  *  1 *MT 
N  =  10* I 

IF (N-MAJ 51*51.50 

50  N~MA 
L^l 

GO  TO  52 

51  L  =  N-9 

52  SSN*0. 

SSNSe0 • 

MM”N-Li-1 

EM-MM 

DO  55  J-L.N 
SSN*SSN>SN ( J ) 

55  SSNS-SSNS+SNI Jl**2 
CHIS«<EM*SSNS/SSN-SSN|*TU ) 

DO  60  J* 1 . 4 

IF<EM-10. 156.57.56 

57  P*0(J) 

GO  TO  56 

56  P*(EM-10.)/10«*<C( I)-B< I) >+B(I) 

58  IF  (P~CHISJ61.60.60 

60  CONTINUE 
J*5 

61  GO  TO(71.70.69.70.71J »J 

70  CONTINUE 
GO  TO  71 

69  PUNCH  200 

PUNCH  6.KEP1.KEP2.KEP3.CHIS.MM 

1-1 

IF  I  JC.K-9)  317.319.317 

71  PUNCH  200 

PUNCH  7.KEP1.K.EP2.KEP3.CHIS  »MM 

1*1 

I  F  ( ICK-*9  1317. 319.317 
317  I F  I  K.K-6  )  102.300.102 

319  PRINT  320 

320  FORMAT ( 16HREAD  NEW  PROBLEM  j 
PAUSE 


GO  TO  302 
215  KK-9 

GC  '0  12 

CO  TO  12 

:  CQRMAT <2XI4tlX 14*1X15. 1 XE10.4 * 18HSAMPLE  REJECTED  0NI5,^H  OBS > 
..•■•>0  f-ORMAT  ( 28MSYSTEM  DAY  TIME  CHISi 

^  ‘0RMATI2X!4*1XI4,1XI5*1XE10.^*18HSAMPLE  ACCEPTED  0NI5*4H  C3Si 
10  4  •-ORMAT  (  21  H«  TIMES*  ARE  NOT  EOUAL  > 

20.;  f  ORM  A  T  (  4  OH  * 

132H  »I 2) 

:«>!  PCRWATl  I<»,  12  »3E9,3  ) 

RCO  ; ORvAT ■l4,t3,!4«I5*!2«2i3,fr8.0*F7.2*F8.0*F?#2*l4,£9.3,F5.2i!2 
f.MD 


tabLe  xxi  output  ch!S6 


<0  Is 

o>  <n 


id  in  uj  o 
Z  m  X  o 
•—00  —  0' 
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Q  —i  Q  r* 


X  ^  X 

UJ  r-«  UJ 


«/) 

>- 

UJ 
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</) 


/pa-  XXI!  Dcryv  I  SOURCE  PROGRAM 


OFCAr  ! 

'  Z 0  Re. AC  100.NYX 
PUNCH  XOO.NYX 

-JlAO  101  *  l  P  *NRl.  M*XI  ♦  XLAM2 .XL AM3 
PUNCH  101  .JP.NRt  M.Xl  *  XI.AM2  ♦  XLAM3 
IP«NRLM~I >340.300.340 
300  READ  100. MYX 
PUNCH  100.NVX 

I F / NYX  +  3 ' 300*301*300  ,  , 

■O'  RfAp  102,lD.LT.;PM.NYX 
PUNCH  102*1.0*:.  T.l  PM.  NYX 
ii  j F :NYX <1 i 10.302.10 
10  READ  100* NYX 
PUNCH  100. NYX 

co  TO  ; i 

OCT  RF.AO  101.  (Al,!A2,IA3.IAi,:AR.!AE,IA7,01  .B2.B3.B4.!BE.0.BE.N*'i 
I f i XL ) 30**500.304 
5  0  0  Y  L  -  D 
3  0  4  I  S  M  ~  I  A  1 
1  GO- I A5 
I G  -  i  A  6 

IF (wD: 20.30: *20 
?Oi  L  D  -  1 A  3 

LT* IA4 

20  t A4* ! A4-«0» i IA4/100 ; -LT+40* ( LT/ 100> 

XUA ] 3 1A3-LC  I 

HUA?= I A4 

25  T  .a XUA  1*1  440  •  -*-XUA  2 
l  T0*B1*F.XP<  XL  *T' 

c  i  INCH  104.  i  Ai  .  JA2  ,LD.IT,  I  A5,  I  A6,  I  A"7*  JB5,ATO,B?  .R?  »B4  *  X!..  ,B*  Nv  i 

PL  AO  105,  !A1 , IA2 . IA3 .IA4.IA5. !A6. IA7.B1.B2 .B3.B4, f B5.0.B5  *NY> 
!  f  <NYX>9',  310.309.310 
3 ! 0  l -  (NYX-6  >312,311.31 2 

2  1  2  :  p  *  ’  SM-  j  A  1  •;  304 , 3 1  3  ♦  30  4 

3  13  ,*  F  {  :  G-  ■'  AC  >304,20*304 
30?  PUNCH  1C0.NYX 

i f <  j  CD- 3 >  r 31 .502.502 
'•O1  PRINT  30* 

304  FORMAT  i  2  "HLOAD  THIS  PROBLEM  WITH  BY,Th> 

50*  PPIH7  350 

?  3  0  f  CRH  AT  i  i  CM..  CAD  NEW  PROBLEM  >  • 
r  r.  ij r-; 

•:-o  to  : :  o 

<  i  1  PUNCH  }  00 . N v  x 
50  T0  10 


340  PRINT  333 

333  FORMAT  1 20HLOA0  DECAY  11  OR  111) 

PAUSE 
GO  TO  320 

104  FORMAT! 1 4. 1 3 , 1 4, I  5 . I  2 .3 l 3 .E9 . 3 ,F7.2 .F8 .0 »F7 .2  .E9. 3 »F5. 2 . 1 2 ) 

105  FORMAT(I4.I3,I4»I5.I2,2I3»F8.2»F?.2*F8.0.F7.2.l4tE9.3,F5.2.12> 

101  FORMAT1I4.I2.3E9.3J 

102  FORMAT1I4.I5.14.59X.I2) 

100  FORMAT ( 40H  • 

132H  .12) 

END 
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TABLE  XXIV  DECAY  II  SOURCE  PROGRAM 


C.  DECAY  11 

DIMENSION  SST ( 3 .30 ) 

50 2  FORMAT i 3F 10.2  I 
i'O  603  1=1,30 

503  READ  502.SSTU »I } .SSTI2.I 1 ,SST<3,I > 

DIMENSION  $N<60> • T(60) ,X1 (60) *X2 (60 ) ,KX(60> 

1  READ  200. NYX 
PUNCH  200. NYX 

READ  219.NP.NI .XLAM1 .XLAM2.XLAM3 
PUNCH  219.NP.NI .XL AMI .XLAM2.XLAM3 
1F*NI-2)217. 202,217 

217  PRINT  218 

218  FORMAK 17HNOT  DECAY  II  DATA) 

PAUSE 

GO  TO  1 

202  READ  200. NYX 
PUNCH  200. NYX 
IF(NYX+3)202, 602.202 
602  READ  203.LD.LT .NP.NYX 
PUNCH  203. LD.LT. NP.NYX 
206  READ  200. NYX 
PUNCH  200. NYX 
IF  (NYXt-1  1206. 207. 206 

20  7  READ  70h.  IA1 . 1 A2  .  I  A3  .  I A4 » I A5  » I A6  .  I A7  .SNBAR  .B2  . YT  ,B4 »  I B5  . D  .66  .NYX 
215  IF (LD)209.208.209 

208  LD -  I A 3 
L  T  *  I A4 

209  1*1 
ISMP-IA1 
IG=IA6 
1X2= IA2 
IX5=IA5 
XX2*B2 
YYT»YT 
X4~B4 
IY5= IBS 
X6  =  B6 

GO  TO  20 

3  READ  704. IA1 . IA2. I  A3, IA4.IA5.IA6, 1 A7 .SNBAR .82  »YT»B4,IB5.D*B6,NYX 
I F ( NYX+  6 1210.16. 210 

210  IF(NYX+9)212»16.212 

212  1F( IA1-ISMP) 16.213.16 

213  IF (  IA6-IG) 16.20.16 

20  IA4*  IA4;-40»  ( IA4/100)-LT+40*(LT/100) 

XUA 1  =  I  A3*-LD 


XUA2“ 1 A4 

25  T ( 11 -XUA1M440.  +  XUA2 
KK  (  I  )r.  IA7 
SN ( I l =$N8AR 
I  °  I  ♦  1 
60  TO  3 
16  N=l«l 

PRINT  1  20.N 
DO  57  1  =  1, N 
XI ( I ) *EXP ( -XLAM1 * T  <  I  )  ) 

57  X2  (  1  ) *EXPI-XLAM2*T  t  I  )  ) 

X1SQ*0. 

X 1 2 -0 • 

X260-0, 

DO  30  I  c  1  »N 
X 1 SQe  X 1 SQ+X 1 (  |  )  *X  2  S  I) 

X2SQ=X2SO+X2  t  !  )  *  X  2  t  I  '• 

30  X12*X1?+X1 • I )*X2<  I  1 

CU  =X2SQ/! X1S0*X2SQ-X12*X12  , 

Cl 2=- ( C 1 1 *X12  )  *X2 SO 
C2  2  =  X  1  SO/ ( X 1 SQ*X2  SO- X 1 2*X 1 2 ) 
C12=-(C22*X12 1 /X1S0 
Z1  =  0. 

22=0, 

00  40  I-l.N 

Z1  *Z 1  +  5N ( I »  * XI ( 1) 

40  Z2=Z2+SN< I ) *X2 < I ) 

A01=C11*Z1+C12*Z2 

A02=C12*Z1+C22*Z2 

SMSO-O, 

DO  50  1*1. N 

SNBR" AO  1 *X 1 ( I  I +A02*X? I  I ) 

50  SMSO*SMSO+(SN( I ) -SN8R ) * ( SN (I J-SNBRJ 
N*N~2 
£N-N 

IF (X6- 1.96 >403.402. 40 3 

402  1=2 

00  TO  410 

403  IFCX6-1.65 1405.404.405 

404  1=1 

00  TO  410 

405  1=3 
410  J=N 

IF  lh-30>407-407,408 
408  J=30 
407  S T  =  S S T  <  1 

XLA01=A01-ST*SOPT  <  Cl  1 *SMSO/EN> 

XUA01  =  A01  +  ST *SOR T  .  C  1  1  *SMSQ«/FN  » 


XLA02*A02-ST*$QRT ( C22*SM5Q/EN  ) 

VUA07.x  A0?  <ST*50RT  (C22*SMS0/EN) 

NY  •*£ 

PUNCH  6 1 3  *NY 

PUNCH  501 »A01 ,XLA01,XUA01 
PUNCH  501«A02*XLA02«XUA02 
Nv«? 

PUNCH  5 1 3  *NY 
N*N  +  2 

DO  400  l  s  1  *N 
E 1 a-XL AMI * T  I ! ) 

E2-=-XLAM2*T(  ]  ? 

ABOX* AO  1 *E  X  P  f  E 1 l+A02*EXP(E2 ) 

RAT  I O-SN ( I ) /ABOX 

AOIC*RAT!0*A01 

AO2C=RA"IO*A02 

M-ISMP 

NY-0 

PUNCH  204, M, IX?,LD,LT,IX5,IG,IOCr I ) , I Y5 ,A01C ,XX2 ,YYT ,X4  , XL AMI  ,X6,N^ 
aOO  PUNCH  204, M, 1X2, ID, IT,IX5,IG, ICKI  I)  ,  I  Y5  »A02C  ,XX2  •  YYT  ,X4  ,XLAM2  ,X6  ,N1 
IF  (NYX--6  »  240,250,240 
2  *  0  IF /NYX+9i215»230»215 

230  PRINT  231 
PUNCH  200.NYX 

231  FORMAT l 13HREAD  NEW  DATA! 

IF ( 1X5-51701,702,702 

70?  PRINT  705 

'705  FORMAT!  14HL0AD  WITH  BICTHl 
701  GO  TO  1 
250  PUNCH  200, NYX 
CO  TO  207 

5  ?  ?  FORMAT |40X,3?X,  12 ) 

501  F C-RMAT  (  3E  12 •  5  ) 

200  FORMAT (40H  ♦ 

1 3?H  1 1 2  I 

?»o  FORMAT! 14, 12  ,3F9,3 I 

703  FORMAT! 14, 15, 14, 4 OX, 19X, 12) 

120  FORMAT  <  3HN  «,I4> 

204  - ORMA T ( I4,I3,I4,15,12,3I3,E9.3,F7,2*F8.0*F7.2,E9.3,F5.2,I2) 

704  FORMAT ( 1 4 , 1 3 , 14 , I  5 , 1 2 ,2  I  3 ,F8 ,2 *F7.2 ,F8 .0 *F7 .2 , 1 4 *E9# 3 ,F5 ,2 , 1 2 ) 

END 


TABLE  XXV  OUTPUT  DECAY  II 
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XXVI  DECAY  III  SOURCE  PROGRAM 


DECAY  III 

DIMENSION  SST  i  30  !  »X._M<  3  )  *XlA<  3  )  *XUA( 3>  ♦  AO  *  ?  ) 
DIMENSION  X3  (  20  !  »  A  I  3  »  3  ) 

D  ? MENS  1  ON  SN(20) *  T  f  2  0  ) *  X 1  (  ?  0  >  »  X  ?  (  2  0  1 
502  r ORMAT  S  F 10*2  t 
DO  503  1-1*30 
M> 3  READ  502*SST<n 
l  RF  AC  200*N'X 
l  -  l 

PUNCH  SOO^N^X 

READ  2l9,NP*Nf  ,*LMU  *  *XLMf 2>  ,XLM{? ) 

PUNCH  219,NP,NI ,XLMl U ,XLMf 2) «XLM(3) 

206  READ  200* NY X 
PUNCH  200 *NYX 
*F«NYX-*1  2  SI  ?  *3  *61  7 
St7  '  F  '  N  v  X  ■*  3  i  2 06*?0?*20 6 
2C2  READ  20:  »!.D*LT*NP*NYX 
PUNCH  203*LD»LT*N°*NYX 
GO  TO  20fc 

216  l?  < LD i  209  *  208  *  209 
'  0  H  i,  D  i  1  A  3 
’*  T  -  1  A  A 
209  ISMF^IAl 
i  "1 

JG-^IAS 

iX2iIA2 

iXS---!A5 

KU!A? 

/  x  2  -  e  2 

YYTi-YT 

xA«e<* 

I Y  5  -  I B  '3 

,C  *0  20 

3  READ  9Q6„IAl«IA2*!A3»IA-*tIA5»IAc!»IA?»5NBAR*B2.  -  YT»6 
IF  t |-\ >560, 215*560 
:<0  IF 'NYx+S>210»16*2l0 
2  i  3  JP'.NrX‘9»212,16>212 
2!2  t  F  i  »  A  1  >-•  I  SMP  ,16*213*16 
■  1  jF  I  IA6-IG1  16*20*16 

-.•>  IA4-IAA-40*  (  IA4/100»-LT*60»K.  T/100» 

U'Ai 1  }  *  |  A  3-S.D 
XUA i 2 i  » L A<* 

f ; I ) -XU A i 1 >  *  1640* >XUA i 2  } 

SNI  i  I  -F.NBAR 


1*1  +  1 
GO  TO  3 
16  N-I-l 

PRINT  120. N 
702  DO  57  1=1, N 

XI  (  I  I  -EXP{-XlM<1 ; ( i j , 

X2 <  I  J *EXP< - XLM  5  2 i *7  III) 

57  X3(I i=EXP<-XiM(3>*T( 1 ) j 
DO  600  1  <=1 ,3 
00  600  J=1.3 
600  AlltJhOi 

DO  30  I l  ,N 

A  < 1 *1 ) *A( 1  ,1 t+Xl < 1 J  *X1 i  I  | 
A;2,2»*'A(2*2)+X2<  I  l  *  X  2  (  1J 
A(3,3)-A(3,3)+X3( I)*X3( 1 » 

A ( 1 »2 ) * A  t 1 » 2 ) +X 1 ( I )  *X  2  ( 1 ) 

A  (1 ,3  »  =AI  1  »3)+Xl(  I  l  * X 3  f  I 1 
30  A(2,3l*A(2,3j+X2< Il»X3( n 
A  (  2 , 1  >  -  A  <  1  •  2  > 

A  (  3 « 1  i  -  A  f  1 » 3  1 
Ai3,2)=Af2«3> 

DO  170  1=1,3 
DUMP  *  A  <  1,1) 

A  <  I  •  I  )  -<  1  • 

IF (DUMP) 130,110,130 
110  PUNCH  131 

131  FORMAT ( 20H0CAN.T  INVERT  MATRIX) 

STOP 

130  DO  140  J*  1  » 3 
140  A ( I , J ) « A ( l , J  > /DUMP 
DO  170  IC  =  1  *  3 
I F  ( K.“  I ) 150,170,150 
150  DUMP=A ( K ,  1 ) 

A?ft.I )-0, 

DO  160  j-1,3 

160  A(K«J)>A(K,J) -DUMP* A  <  I  »  J  ) 

170  CONT inue 

Zl»0. 

22*0, 

Z3»0. 

DO  40  I = 1 *N 

zi*zi-»sn<  n *xi  m 

22  *Z  2  *  SN (  I  l*X2<  I  ) 

40  Z3  =  Z3fSN(  I  J*X3< I ) 

A0fl|eAil«l»»Zl+A<2'll*Z2*A(3,l)»23 
A0'2l=All,2)»Zl+A)2*2MZ2*Af3,?n)Z3 
A0»  3  *  =  A<  1  »3>«Z1-A  )  2,3  >  *Z2  *-A'  3,3  )  *Z3 
SMSQ-O, 


E'  f;o  >0  :-i,n 

W**-:  o*  \ '  ** 1  ?!  >  *  AO;  2 \>  x? '  I » -An.;  2  j  */3»  j  • 

Rc.Lf.fi  •  ABS(  GNU  »  l /SNOR- „? 

U  iRELfh:&0#-i0*  *00 

GO  SMSO  -5MS0  '  •  !  .f*N  •  J  *  -  SNf-o .  *<  MU'-  SM6R  V  j 

v//  to  M  N  *•  3  , 

1  1  :  N  -  N 

■•‘10  ,!  -to 

;•  fN~30<«CU»«0?»*08 

o  p  ;•  -•  3  o 

ii;  '-o >  .s t •■» s r> T /  ji 

CO  571  Jsl*3 

A  <  J ,  ,1 }  vSt  *S0R1  t  A ;  J  «  J I « SMOO/EN I 
XL  a  :  J',  v  t.C. c 

•  j  !  r(:  a  :  j  , AO »  >  i  +  A  ?  v  ♦  J  i 

?4  V-r> 

p(jtovn  &  i  ? ,  to  y 

00  0/2  J-;,1 

:  •  ?.  PUNCH  00J  »A0w?4*LA(  J5*XUA*  J5 
f  1  v «  7 

PUNCH  r»:3,NY 
N-N  •? 

!>0  uof)  T  '  1  ♦  N 

ABOX^AO*.  1  f  -VI  (  I  }  «  A0(  ?  )  *X2  (  1  !••  AOi  si  <  y.M  t  1 
RA r lO-SNi » > .' ABO* 

MUSMP  * 

N  V  --  0 

f-0  400  .>=1,3 

A  <  *  J  f  »  R  A  T  !  0*  A  0  4  J  »  ' 

-  0  3  PUNCH  ?;•<->«♦  1X2  »Lf'  ,1.  r  >  »  >.3  ,  fG,VK  ♦  T  Y0,A  i  J*.)>  ,  XX2..YYT  ,A4  ♦XLM{  J?  ,Xf =  ,* 
If JtoYX+61 2a0,200*;4O 
2  4  0  tr  9 , 215,2  30*215 

730  PRINT  231 

PUNCH  200»NVX 

231  FORMAT ( 1 3HREAP  NEW  DATA» 

Jr  5  !X5<-0)901. 902,902 

!07  PR  ?  NT  *>05 

•  jQc;  FORMAT  U^HICAD  W 1  »  H  8  K  T  h  } 

>31  GO  T0  i 
TOO  PUNCH  iOOfNYX 
1  \ 

;.o  TO  3 

00  Hi  l  J*  *  ,Nr 
S.N(  JJ  -Gfl*  i‘  ~l  i 
TOT  T  i  -T  «  Jf  1 1 
N.  N  -1 

■::  co  to  ?:*? 


513  FORMAT (40X.32X.12 ) 

501  FORMAT (  3E 12 .5  ) 

200  FORMAT (40H  • 

132M  *12) 

219  FORMAT! 14,12, 3E9. 3) 

203  F0RMATH4, 15,14, 40X,19X»I2) 

120  FORMAT ( 3HN  a, 14) 

204  FORMAT! 1 4 , I  3 , 1 4 , I  5 ♦ 1 2 ,3 1  3 ,E9. 3 .F7.2 ,F8 .0 ,F? .2 ,E9.3 ,F5.2 »  12) 

904  FGRMATU4,I3,I4,15,12,2I3,F8.~,F7,2,F8.0.F?.2,I4*E9.3*F5.2,I2) 

'  END  .  ■  ' 
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TABLE  XXVIII  VAR  I  SOURCE  PROGRAM 


C  VAR  I 

DIMENSION  T(3,31) *X(200> 

DO  1  1=1,31 

1  READ  103,T<1.1 ).T<2.1 ) »T<3»I> 

103  FORMAT ( 3E10.5 ) 

10  READ  100. NYX 

IF (NYX+1 J 10*11*10 

11  READ  101.  IA1,IA2,IA3,IA4,IA5»IA6,IA7,IB5.B1.B2»B3.B4,XLL,B6,NY> 

12  N  =  1 

I SMP  = I A 1 
JGD= IAS 
IGa I A6 
X ( N ) =B1 
I PROB= I B5 
I SY= I A2 
XEP=B6 

17  READ  101.  IA1.IA2.IA3.IA4.IA5.IA6.IA7.IB5.B1.B2.B3.B4.XLL.B6.NY) 
1F(NYX+6J 13,50,13 

13  IF(NYX+9) 14,50,14 

14  1 F ( IA6-IG)50, 15.50 

15  IF( IA1-ISMPJ50.16.50 

16  N=N+ 1 

X ( N ) =Bl 
GO  TO  17 
50  M=N-1 
XBAR=0. 

SXSQ-O# 

DO  18  1  =  1, N 
XBAR  =  XBAR  +  XU) 

la  sxsossxsQ+(xm*x(  m 

XN=N 

SSX=XBAR#XBAR/XN 
XBAR  =  XBAR/XN 
I F { N-l ) 20.70,20 

20  S 1GWG= ( SXSQ-SSX ) / ( XN-1 • ) 

SIGWG  =  SORT  CSIGWG) 

SXBR=S IGWG/SQRT (XN ) 

IF (M-31 172,72.21 

21  M=31 

72  IFIXEP-1. 65)63, 62.63 
62  1=1 

GO  TO  22 

63  IF(XEP-1.96)65,64,65 

64  1=2 

GO  TO  22 


65  1*3 

22  XL-XBAR-SX8R#TU ,M) 

U-XBAR+SXBR*Tl I »M) 

PCERR*><SXBR«T<  I *M ) /XBAR 1*200* 

Nl* 1Q0*N-N 

N2*100*N/4-N 

N3«100»N/9~N 

N4*1Q0*N/16-N 

N5*l Q0*N/25-N 

N6*l00*N/36*-N 

N7*l 00*N/49“N 

N8*100#N/64-N 

N9*100«N/8l-N 

N10  =  0 

PUNCH  110 

PUNCH  105.1 PROB . I SY  .  1 SMP ♦ 16  *N 
PUNCH  111 

PUNCH  106.XBAR.SI 6WG .SXBR.PCERR 
PUNCH  112 
PUNCH  107  .XL  »U 
PUNCH  113 

PUNCH  108.N1  .N2.N3.N4 .N5 .N6 .N7 .N8 *N9 *N 1 0 
IF (NYX+6 )61»11.61 
61  lF(NYX+9)507.80.507 
61  IF (NYX+9 ) 12*80*12 
70  PUNCH  110 

PUNCH  105.IPROB.I5Y, ISMP, !6,N 
80  PRINT  200 
PAUSE 
60  TO  10 

507  IF  <  JGO-8 1509*508*509 

508  PRINT  510 

510  FORMAT  I 15HLOAD  WITH  CALlB) 

60  TO  *2 

509  IF( JGD~9ll2.512.12 

512  PRINT  513 

513  FORMAT  I 14HL0AD  WITH  fLUX) 

60  TO  12 

200  FORMAT  I 16HREAD  NEW  PROBLEM) 

110  FORMAT  (  3X  *  7HPROBL  EM  ♦  4x,fcHSYSTEM»4X6HSAMPlE,2X  »  8HGFOME  TR  V  * 

16HNUMBER I 

107  FORMAT ( 2F 10* 1  I 

108  FORMAT ( 10161 

111  FORMAT < 10HMEAN  COUNT , 7X . 3HS 1 G ,6X .4HSXBR ,2X , 1 1 HPERCENT  ERR) 

105  FORMAT  (  4X  *  1 6.  »4X  *  1 6  *4X  *  I  6  »4X  »  1 6  *4X  •  I  6  ) 

106  FORMAT{F10.0.F10.5*2F 10.4 j 

101  FORMAT  f  I4*I3*I4*I5»!2*3!3*E9.3*F?,2»F8.0«F7.2*F.9.3*F5«2»  12) 

112  FORMAT ( 2 OH  LOW  LIMIT  u»  LIMIT) 


1 00  FORMA  r i 

132H  '  .  4 ;  *  v 

113  ?  OR  M  A  f  <  3  x  ,  3  h  •  I L  *  3  x. »  3  M  •  2 1  >  3  a  4  3 1 1 .  ‘  l  4 3 3m**E ♦ 3* 4 3m* if , 3  a  .3ii.SC«3i 
•  I  »N«  n. « 3*  *  3H*«t  *r<  « ?»I4  9C  *3X3‘tl  ♦£. ; 

END 


UNO*  4  ’  3  *  V U  *  S •  GM*v  *  V  t‘ > *  •  » S  !  S 
4  3  3  *  ORMA  T  (  3£9*3*$M  \  A. 

c:::-,;-  ru  ••.:>*  ; u  .  '  : 

» 3<*  *  Of'MA  7  (  » 6H3AMPL6  S  IN  GR'V'P  * 

PUNCH  73  3  « l  !  X  )  «  i  <  ?  I  •  1  J  3  !  »  I  i4  »>  M  5  J.*.!  ••  6  )  *  l  <  7)  »  I  (  6  »  « I  >  9  l  ♦  I  ■  X  0  ) 
.735  -QRMAT  1  1CIM  •  >';■ 

PUNCH  407  ,  ,ak:  :  V ■■•  :V"- 

407  fORMAfi  1 3MPERCEN7  ERROR’ 

c.jN'H  1 1 3  * 2 Z  M  !  *  2  2  t  ?  '  *  2  2 «  3  '  ♦  Z  2  ’  *-  »  » 2  Z  ;  E  i  « 2  2  if  6  ?  *  2  2  i  ?  5  ♦  2  Z  »  8  * »  2  Z  ( ‘ 
!  22  {  10  »  .  ,-0 ;  . 

¥:0t  PUNCH  609 

VO 9  FORMAT; ifHMR  SMf*L  $  NEEDED. 

S-P# 

OUNCM  113  »xm;  f  VM2  ,XM3  *XM4  »XM7  .XMt*>M?  fXM8»XM9»G 
XF  fNYX  +  t  ^  700 .300*700  .-..A-A 

?00  CONTINUE  ^vv‘:V 

?R  TNT  PEE 

>*«*  FORMA  T  (  16H.CAC  NE #f  PROBLEM  .* 

PAUSE  ’  7:  .'V. 

GO  70  H. 

104  FORMAT,  i  4*  >3*16*1  5  *  ?  X  « J  7 «  3*  *  ;  3  *h'9*  3  »  22  X  *L9*  7  *5X  ♦  T  3  ‘ 

200  FORMA  T  j  40H  • 

1 32H  *I2l 

?01  -09MAT  2i  *>.40X*?2X.I2  ; 

JOG  "  ORM  A  T (  ;  4HR  IS  IMAGINARY* 

13  FORMAT  *  lCr -’.l  > 

END 
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TAR'.t  xr>,  PRE  f  /AR  ;  !  ;iOuR'.:t  PROGRAM 


C  F»EP  r op  va«  It 

v  M3X  25  :->.OME  [RTFS 

DIMENSION  tCll9i  *EC2<9i  *EC3'.9>  .  IGl2S  »  *STt9 
Li  0  SOI  .  -•  n  9 

SOI  READ  2lQ*CCl'Ul ♦EC2. JJ,EC3(J> 

DO  302  »=1 ,9 

302  READ  303»ST(Jj 
203  ° CRMA T  ;  r  1 C  j •»  '* 

I  K’~i 

DO  SOO  !  =  1  *  2  *> 

900  i  G  »  M  ^ 0 
NGOS  -  0 

10  Rt*r'  201  ♦  N v  j> 

PUNCH  r01*NvX 
READ  200,  r  PROP  »NRi.  AM 
PUNCH  OGOtIPROB.NPLAM 
READ  ?0?»NGC 
PUNCH  202  «.NGC 
?2\  READ  201*NYA 
PUNO  2  0 1  ,  N  Y  X 
• FtNYX-7,221 ,21*221 
?i  READ  20“ * 1A6,NY< 

i  f  *  r«vx— 6 «  aoo*'*oj  » coo 

*CI  READ  2CA«JA6,Nry 

IF  <NYX-?»401*21*‘r0i 
GOO  IF  f  N  y  X  ♦  6  <  f-  0  ♦  3 1  *  r  0 
81  PUNCH  204*IAo*NyX 
GO  TO  21 

50  IT  I N  y  x  -f  9  I  30*  60 1  *30 

30  SO  T 0 1 3! *32 ,33» »KZ 

31  PUNCH  £04  *  I A6 »NYX 
READ  2Q4,IA6*NYX 

GO  TOI *60, 700*61 ) vNRI AM 

51  RE. AD  20*,iA5,NYX 
GO  TO  700 

32  READ  2Q4.1A6.NYX 
PUNCH  204  * IA6.NYX 

GO  !Ot 7c0.21 *61 > *NRLAM 
62  READ  2 0 *  » I A  6  *  N  Y  X 
GO  70  21 

33  READ  204, J  A  €  *  N  Y  X 
S  3  RF  AC  ZOi-,  IA6  «NYX 

PUNCH  204,  I A  £  »  N v X 
00  To  21 
SOD  KK< l 


00  701  I >1 *25 
IF  <  IA6-IG( I ) ) 701  *  702  #701 
702  KK-0 
701  CONTINUE 

IF<KIC)704. 21.704 

704  00  706  1*1.25 

I F ( IG  < I )) 706.705.706 
706  CONTINUE 

705  1GIIIMA6 
NG0S*NGQS+1 
GO  TO  21 

601  PUNCH  204.IA6.NYX 
PUNCH  208.NGOS 

208  FORMAT ( I  3 ) 

IDF«NGC-1 

DF-IDF 

PUNCH  209 .OF .ECl ( IDF) ,EC2(IDF) .EC3( IDF) 

00  503  J-l.25.5 

I*J+1 

X«Jf2 

L-J  +  3 

M»J+4 

503  PUNCH  504.IG(J).IG(I)»1G()C).IG(L).IG(M) 

504  FORMAT <51 5) 

PUNCH  303 .ST  < IDF ) 

IF ( KZ-NRLAM ) 60 2.760. 760 

602  ICZaXZ-t‘l 
PRINT  780 

780  FORMAT ( 11HREREAD  DATA) 

GO  TO  70 
760  PRINT  603 

603  FORMAT ( 16HREAD  NEW  PROBLEM) 

PAUSE 

GO  TO  1 

209  FORMAT ( 4F 10 .6 ) 

20*  FORMA  T ( 18H  .I3.19H 

1 32H  .12) 

201  FORMAT  <  *0H 

132H  .12) 

200  FORMAT  < 14. 12  »34H 

1 32H  .12) 

202  FORMAT < I 5.35H 
132H 

210  FORMAT ( 3F 10.3 ) 

END 
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I  V  W\U  VAk  11  SOURCE"  PROGRAM 


C  VAR  II 

DIMENSION  IGl25),FV19)*I<10) *SUM< 10) *SUM$Q< 10)»1<10>*ZZ(10) 
DIMENSION  FT(102) .XNBAR(IO) 

41  READ  501.NGOS 

READ  502,DF,EC1,EC2.EC3 

501  FORMAT (I  3) 

502  FORMAT  1 4F 1 0 • 4 ) 

DO  503  J®  1 *25*5 
1  I  » J+  1 

K-J  +  2 
LL« J+3 
M*  J  +  4 

503  READ  504,161 J) ,IG< II) ,IGiK) ,IG(LL> *IG(M) 

504  FORMA  T i 5  I  5 ) 

READ  303, ST 

303  FORMAT1F10.4) 

IDF-DF 

DO  910  J-1,102 

READ  2 02  » p V  (  1 J ,FV(2) *FV(3) »FV(4) »FV(5; * FV ( 6 ) *FV(7) ,FV(8) *FV19) 
910  FT  C JJ-FV1  IDF) 

202  FORMAT ( 9F7#2 ) 

1  DO  700  KX-l.NGOS 
I  F  ( XX- 1)827,817,827 
827  PRINT  254 

254  FORMAT (11HREREAD  DATA) 

817  READ  200, NYX 

READ  940, IPROB.NLAM 
940  FORMAT ( 1 4 , 1 2 ) 

READ  201 ,NGC , I PROB  »N YX 
XM-NGC 

IF(NGC)211, 212*211 

212  PRINT  213 

213  FORMAT (10HVAR  I  DATA) 

PAUSE 

GO  TO  1 

211  DO  215  J* 1  , NGC 

215  READ  201 , I ( J) , IPROB 

216  READ  200  »NYX 

IF (NYX-7)216*  300 ,216 
300  DO  230  «J«1»1Q 
LIJJ-0 
iUMi j)-0, 

230  SUMSOl J ) *0, 

S»Q. 

SSQ-0. 


ev  ?>-»?>• 

FT5*0. 

+  •{)* 

YM'Jl-Oi 

«‘-0* 

S  I G  i  =»  0  « 

23  S  Rf  AO  104  ♦ I A 1 ♦ l A2  i I  A3  » I A4  *  FAStlBStBl »XLAM»N: X 
ir(NYX+9i225»219#225 
22S  IT *NYX*6'.226,219,?2s 
22£  DO  232  J= 1 « NGC 

I f J I A  1  - F i J> )232»233«232 
23  3  I?  (  !A6~!C(KX.n23f>t?24  »2  3S 
232  CONTINUE 
224  sum: J j -SUM{ 3 J+Bl 

SUM3Q  *  J  l  -  SUM  SO  (  J  |  -*-B  1  *  8  1 

L  l  J  i  B  L - J )+  1 

iDYwr  f  A3 

ITEM- 1 A4 

IS?S»IA2 

?i  AM*  Xt  AM 

GO  TO  2  3C’ 

21?  DO  241  n=\»NGC 
S*S*SUMiK) 

5S0-SSQ  SUMSOtKi 
*NI =L (XI 
T - T  +XN 1 

SIG1*S IG1+1 */XNl 
SQS«<SUM(K>*SUM<K  < }/XNj 
EVS«EV5>M5!jM5Q(K|~SQS  ) 

AN?  At* «  r  j  ,  ,  'YK<  f  *  SUM  t  <  i 
yMUI*YMU1  +  XNBAR.'K  » 

241  YMU-»t  U/XM,*YMU1 
DO  16  X*  1  »NGC 

l<?  R2rR2*  (  XNBAR  <X  j  -  YMU  t  *  <  XNBAR  l  K  *  -YMU  I 
OS*S#  T 
ETS-SSQ-C 
EM$  »E  TS*-  E  VS 
ETD*T-l. 

FHD**AM-1. 

EVD^ETD- EHD 
EHM-EHS/FHD 
f VM-EVS/EVD 
ETMeETS/ETD 
EE»FHM/£VM 

R1 -R2.XMM 
I E : R 1  i  11  1 1 2  *  1? 

11  PRINT  1  OS 


GO  TO  700 
12  RiSQRTlRl) 
lobs-evd 

IF (LOBS* 100  >902.902.903 

902  TC»FT(LOBS) 

GO  TO  1000 

903  1 F ( LOBS- 4 00 ! 904  *90 A ,905 

904  TOLOBS-lOO 
TC  »T  C*£C 1 
TC-TOFTUOO) 

GO  TO  1000 

905  JF(LOBS-I000)906«90t.907 

906  TC2LOBS-400 
TCrTC*EC2 
TC-TC+FTI 101 ) 

GO  TO  1000 

907  TCaEC3 

1000  IF !FP-TC> 100*101* 101 

101  PUNCH  340 

340  FORMAT l 8HVAR  SMPS > 

100  TC-ST 

YEST-RtTC 

Z«100.MYFST/YMU) 

PQ*Q, 

DO  271  X-1.10 
PQ-PQ+,1 

271  2ZnC)-PQ*2 

XM1*XM/0.01~XM 
XM2*XM/0,04~XM 
XM3-XM/G ,  09-XM 
XM4=XM/0,16-XM 
XM5=XM/0 • 25" XM 
XM6*XM/0,36-XM 
XM7-XM/G.49-XM 
XM8=XM/0.64~XM 
XM9=  XM/O • 8  1  -  XM 
SIG2-<EVM/!XM*XM> )#SIG1 
SIG3«XMHR*R-SIG2  I 
IF  *  S I G3  >  14,15.15 

14  SIGMA-O, 

GO  TO  49 

15  SIGMA-SORU  SIG3I 

49  PUNCH  109 

109  f ORMA  T ( 5H  PROB.5H  DA  Y  »  5H  TIME.5H  GE0.7H  LAMDA » 
PUNCH  110, 1 PROB, IDY, I T E M » I G ( XX  I *ZL AM 

110  FORMAT  <  4  l  5  »E9,3 ) 

PUNCH  108 

108  FORMAT (7X»  2HMU ,4X  »  5HS I GMA »5X»4HyEST , IX, 6H SYSTEM) 


P  ONO*  '* .?  3  *  W'J  «  S  '  G  M  A  ♦  v  5  S  T  *  »  S  '  3 
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<•  OP.M  A  “ 

< 3tV*3*iM 

Fir  -  ,;m 

‘  *- 

»  "itt 

-  OiVMA  T 

(  '  6  *■*  o  A  M  P  |_  t  S 

IN  C-RO>‘P, 

t.jnch 

T30  « I  !  I  1  *  !  i 

2 » *  I  :  3  :  * ! ;  <4  > .  i  >: 

*(‘‘6>*I<  7  )  *  I  (  &  > 

I  i  9  >  ♦  I  ’ 

10) 

7  35 

;:OkMAT 

UCii! 

-U'vCH 

40  7 

u  0  ■; 

* ORMA T 

i  1 3M?tRC F  NT 

ERROR ' 

-'JN'H 

i  1 3  .  Z  Z  f  1  :  *  2 

Zt?i*22»3>*ZZ'~ 

>  ,21 ■ 5  :  «ZZ < 5 • *2Z * 

7  5*2 2  i  e ■ 

■  *  Z  2  ( 

Pi.'NCH  *09 

VO  7  POHPAT;  if-HNR  SMPL  0  NS  FDf.O  * 

S  -0. 

°UNOt  113  *XM ;  ,  XM2  *XM3  «XM4  ,XMr  ,  XMt  «  >M  7  ♦  XM8  ♦  XM9  »G 
IF  )  NYX  +  i  >  ’00*300*  700 

700  CONTINUE 
►  •  r t  n  T 

2tb  FORMAT  <16H_C-AC  Nl**  °90Bl f  M  : 
e  A 

00  T0  ti 

)P<*  -ORMft  T  .  •  4*  v  3  ♦  r<»  *T  5*2  X  •  T  ?  «  3>  ♦  :■  3  . 1 9,  3  *  22  X  *L9  *  3  *  5  X  ♦  i  7  *. 
200  ►ORMATj/vOh 

1 32H  *12} 

701  FORMAT  2t*«*CX«?2X*I2i 
!Or>  "7;KMiT(j4HP  IS  IMAGINARY* 

>  1  3  FORMAT  I  J  Cr  -’.l  i 

END 
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ABLE  X  X >  ;  v  BAX  7 H  SOURCE  PROGRAM 


THICKNESS  OF  BACKING  MATERIAL  DETERMINATION 
3  READ  100  »€A  ♦  J GARD 
100  F ORMAT  <F10»3*14> 

I  FUGARD-999  {300.304,300 

300  PRINT  301 

301  FORMAT  f  19HMAX  BE'»‘  A  CARD  CHECK) 

PAUSE 

GO  TO  3 

30*  READ  llO.Nvx 

I F *NVX-1 i 30<*»305,304 

305  READ  102,  JSMP,!SYS,IA5»IG,!PR0B»B1»NYX 
£02  i SM$  = I SMP 
I  «0 
J*Q 
AA-0 
A0*0 

400  ix  UA5^5)404»A03t-40A 
<*03  A0^B1+A0 
1*1  +  1 
GO  TO  500 

404  [X ( JA5-6  1406*405*406 

405  AA«R  1  + AA 
J  »  J 1 1 

GO  TO  500 

40%  IF (NYX+9 >608*501 .608 

608  IF <NYX+6)609, 501*609 

609  PRINT  407 

40  7  FORMAT < 13HNOT  BKT H  DATA) 

PAUSE 
GO  TO  3 

500  READ  102 » JSMP, ISYS, JA5, IG,J PROS, Bl ,NYX 
JF ( I  SMS- I SMP i 501  ,400,501 

501  X‘  I 
AO*AO/X 

r.  J 

aa-aa/x 

7  ] F i j  G“ 1 0  I  10, 11*11 

10  ALFHA=33./<EA**1.33> 

GO  TO  12 

11  I F I EA-l ,45  »  13*14,14 

13  ALPHA*6Q,/ I EA**1 *  54) 

GO  TO  12 

14  AtPHA-200*/ i EA»*6 ,25  * 

12  -LOG i AO) /ALPHA-LOG ( AA ) /ALPHA 

PUNCH  7  05,  I  SM.S  ,  X  F 


IF  <NVX«9l601  *901  #601 
601  |F  tNYx.6)602  #305.602 

901  PRINT  902 

902  FORMAT t 16HREA0  NEW  PROBLEM  > 

PAUSE 

GO  TO  3 

105  FORMAT  < * 1 5HSAMPLE  NUMBER  » , I  4 « 3X ♦ 1 1HTHICKNESS  =.F9.6) 
102  FORMAT ( t 4 ♦ I  3 . 9X • I  2 # 1 3 # 3 X . I  3 #E9 • 3  ♦  36 X  •  1 2  1 
HO  FORMATION  ♦ 

132H  .12) 

END 
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table  XX*V  OUTPUT  BAK1H 
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